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Abstract 

Fiber-reinforced  silicon  carbide  (SiC)  composite  materials  are  important  for  many 
applications  due  to  their  high  temperature  strength,  excellent  thennal  shock  and  impact 
resistance,  high  hardness,  and  good  chemical  stability.  The  microstructure  and  phase 
composition  of  SiC  composites  can  be  tailored  by  fiber  surface  modification,  the  process 
parameters,  and/or  fiber  preform  architecture.  One  process  by  which  SiC  composites  can 
be  produced  is  chemical  vapor  deposition  (CVD).  This  thesis  primarily  focuses  on  mass 
transport  by  gas-phase  flow  and  diffusion,  chemical  reaction  in  gas  phase  and  on  solid 
surfaces,  and  thin  film  fonnation  on  curved  surfaces,  which  are  fundamental  to  the  CVD 
process.  We  highlighted  process  parameters  that  can  potentially  affect  the  structures  and 
properties  of  the  CMCs  using  simple  model  material  systems.  We  also  analyzed  the  use 
of  a  finite  element  modeling  tool,  COMSOL  Multiphysics,  to  build  the  series  of  models. 
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FINITE  ELEMENT  ANALYSIS  MODELING  OF  CHEMICAL  VAPOR 


DEPOSITION  OF  SILICON  CARBIDE 

I.  Introduction 

Overview 

Ceramics  and  ceramic  composites  are  versatile  materials  used  in  almost  every 
industry.  Silicon  carbide  (SiC)  in  particular  is  capable  of  high  hardness,  good  chemical 
stability  across  a  range  of  temperatures,  low  thennal  expansion,  and  very  good  thermal 
shock  resistance.  These  properties  are  made  even  better  by  the  excellent  impact 
resistance  that  can  be  achieved  with  the  addition  of  ceramic  fibers  reinforcing  the  SiC 
matrix  into  a  ceramic  matrix  composite  (CMC).  The  abundance  of  methods  for  achieving 
SiC  production  also  make  it  a  versatile  material.  Silicon  carbide  is  often  produced  from 
polymer  precursors  which  can  be  utilized  in  any  phase  of  matter.  For  example,  slurried 
solid  particulates  can  be  layered  with  ceramic  fibers  or  fabrics  and  then  pyrolyzed;  liquid 
polymers  can  be  infused  via  capillary  action  or  injection  molding  and  also  pyrolyzed; 
gaseous  polymers  can  deposit  directly  upon  a  prefonn  surface.  These  methods,  though 
completely  viable,  are  often  expensive  in  material  costs,  energy  expended,  and  time 
required  [1-3]. 

Such  factors  are  significantly  affected  by  the  particular  manufacturing  process 
that  is  used.  In  this  thesis  the  focus  is  on  the  process  of  chemical  vapor  deposition 
(CVD).  This  is  a  process  that  uses  the  reactions  of  gaseous  precursor  chemicals  to 
deposit  the  desired  material  (i.e.  SiC)  on  a  substrate  or  prefonn  [4].  Optimizing 
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manufacturing  processes  by  experiment  requires  a  great  deal  of  effort  and  cost  to  generate 
more  data  to  fill  the  gaps  in  our  collective  database  o  f  knowledge.  Modeling  and 
simulation  plays  a  strong  role  in  searching  for  effective  solutions  at  lower  cost.  One 
powerful  method  of  modeling  a  system  is  by  means  of  finite  element  analysis  where  the 
large  domain  is  divided  into  a  mesh  of  small  subdomains. 

Purpose  of  Research 

A  goal  of  this  thesis  is  to  gain  an  understanding  of  ceramic  matrix  composite 
processing  by  way  of  computer  modeling.  There  can  be  great  difficulty  during  the 
manufacturing  process  of  CMCs  in  determining  all  that  is  occurring  at  any  given 
moment.  In  situ  measurements  of  material  properties  are  problematic  if  not  impossible  to 
obtain.  It  is  nearly  impossible,  for  example,  to  detennine  the  deposition  rate  at  an 
specific  interior  point  using  traditional  methods  during  the  deposition  process.  However, 
computer  models  contain  all  material  properties  at  all  times  and  locations;  digital 
“probes”  can  thus  be  inserted  into  the  model  and  material  properties  can  be  calculated. 

The  understanding  developed  from  selected  model  material  systems  can  also  aid 
in  guiding  ceramic  design.  Model  parameters  can  be  parameterized  over  a  range  of 
values.  The  resulting  output  can  lead  to  understanding  of  how,  for  example,  temperature 
affects  the  properties  of  the  final  CMC.  Input  parameters  can  then  be  optimized  for  a 
desired  result  in  final  properties.  All  this  is  not  only  faster  than  performing  physical 
experiments,  but  much  cheaper  in  material  and  energy  costs. 
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While  computer  modeling  and  simulation  is  a  relatively  low-cost  and  versatile 
approach  towards  better  understanding  of  a  CVD  system,  it  can  still  be  quite  difficult  to 
generate  an  effective  model.  CVD  processes  involve  macro-scale  material  transport, 
complex  chemical  reactions,  and  sub-micro-scale  kinetics  on  surfaces  that  result  in 
changes  in  material  geometry.  Of  the  various  literature  that  were  reviewed,  two  main 
points  of  modeling  research  focus  were  presented.  Some  authors  discussed  the 
macroscopic  chemical  interactions  and  large-scale  deposition  growth  rates  of  the  SiC, 
many  of  which  uses  a  finite  element  method  model  [5-16].  Other  researchers  discussed 
the  microscopic  changes  in  geometry  resulting  from  the  growth  without  discussing  in 
detail  the  chemical  interactions  leading  to  the  growth.  They  modeled  the  growth  by  using 
various  methods,  a  random  walk  algorithm  being  popular  [17-24].  Two  research  papers 
were  found  that  incorporated  both  aspects  of  this  ceramic  production  [25,  26].  Both  used 
two  separate  modeling  systems,  one  for  each  spatial  scale,  that  were  connected,  but 
neither  used  a  single  unified  model.  In  both  cases,  the  models  would  be  run  iteratively 
and  alternately.  Geometry  changes  affect  gas  flow,  which  can  affect  chemical  reactions, 
which  in  turn  affect  further  changes  in  geometry,  all  at  the  same  time.  As  they  are  so 
interconnected,  a  simulation  strategy  that  can  include,  simultaneously,  both  detailed 
chemical  reactions  and  geometry  changes  is  tremendously  helpful  in  evaluating  any 
manufacturing  process. 

A  secondary  purpose  of  this  research  project  is  to  provide  an  investigation  into 
one  viable  option  for  a  complete  chemical  and  geometric  simulation  of  SiC  deposition. 
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In  searching  for  an  appropriate  modeling  and  simulation  tool,  desirable  qualities  are  that 
it  utilizes  first  principles  to  simulate  multi-physics  problems  and  that  it  be  user-friendly, 
having  preprogrammed  physics  modules  with  automatic  integration.  Due  to  its  module- 
based,  customizable,  multi-physics  architecture,  its  utilization  of  finite  element  methods, 
and  its  general  user-friendliness,  COMSOL  Multiphysics®,  produced  by  COMSOL,  Inc., 
was  chosen  as  the  software  platform.  This  project  seeks  to  explore  the  feasibility  of  using 
this  software  to  model  and  simulate  the  chemical  and  geometric  growth  of  SiC  CMCs  via 
CVD. 

Research  Objectives 

This  thesis  is  focused  on  producing  a  set  of  relatively  simple  models  of  the 
production  of  fiber  reinforced  silicon  carbide  (SiC/substrate).  These  models  are  meant  to 
build  understanding  of  the  chemical  vapor  deposition  (CVD)  process  with  which  the  SiC/ 
substrate  is  produced.  These  models  are  also  meant  to  show  that  it  is  possible  to  build  a 
single  model  that  combines  simultaneously  running  models  of  macroscopic  flow  and 
chemistry  with  microscopic  geometry  changes.  The  specific  SiC  CVD  process  to  be 
modeled  uses  silane  and  propane  as  the  infiltration  chemicals  carried  by  gaseous 
hydrogen.  In  this  case,  the  hydrogen  acts  as  both  a  carrier  and  a  reactant. 

The  process  of  producing  a  model  for  any  complex  system  can  be  made  much 
easier  by  breaking  down  the  system  into  separate,  relatively  simple  parts.  Once  the 
smaller  models  are  generated  and  their  properties  understood  individually,  they  can  be 
brought  together  in  a  unified  whole.  These  sub-models  include  a  model  of  laminar  flow 
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of  gaseous  materials  through  and  around  the  fabric  preform;  a  model  of  the  diffusion  of 
chemical  species,  heat  transport,  and  chemical  reactions  involved  with  the  deposition 
process;  and  a  growth  process  implemented  by  a  deformation  of  the  finite  element  mesh. 
Each  of  these  models  have  its  own  complexity  in  coupling  several  physics  phenomena 
coherently.  These  complexities  are  to  be  explored  by  investigating  the  relationship 
between  initial  variables  and  predicted  results.  How  does  the  input  pressure  effect  fluid 
flow  characteristics?  How  does  temperature  effect  deposition  rates?  Once  these  and 
other  questions  are  answered,  the  three  sub-models  will  be  combined  into  a  overarching 
model  for  the  entire  process.  In  the  event  that  a  working  combined  model  is  achieved, 
higher  levels  of  complexity  and  realism  will  be  added. 

Assumptions/Limitations 

Ceramic  production  methods  are  so  complex  that  it  is  difficult  to  achieve  a  single 
modeling  module  that  will  contend  with  them  all.  Therefore,  for  the  purpose  of  this 
thesis,  only  a  single  process  is  investigated,  specifically  a  chemical  vapor  deposition 
process.  The  models  to  be  used  are  also  fairly  simplistic.  The  substrate  fibers  that  make 
up  a  fabric  preform  are,  for  example,  modeled  as  perfectly  shaped  and  completely  solid. 
True  fibers  have  imperfections  and  porous  behavior.  The  geometry  of  the  initial 
experimental  setups  only  include  the  areas  around  the  substrates  and  ignore  the  inflow 
and  outflow  piping  that  subtract  from  perfectly  laminar  flow  in  the  systems. 
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II.  Literature  Review 


Chemical  vapor  deposition  is  a  process  by  which  matter  is  is  deposited  by 
chemical  reaction  of  gas-phase  precursors.  Chemical  reactions,  often  endothennic,  occur 
in  the  gas  phase  as  well  as  on  the  surfaces  in  the  reactor.  This  can  result  in  deposition  on 
surfaces  of  the  reactor  itself  rather  than  just  on  the  substrate.  Product  gasses  are  vacated 
from  the  reactor  [27].  The  deposition  rate  and  unifonnity  of  a  given  CVD  process  are 
governed  by  the  rate  of  mass  transfer  from  gas-phase  to  substrate  surface  and  by  the  rate 
of  surface  reaction.  At  atmospheric  pressures  with  horizontal  substrate  surfaces,  these 
two  rates  usually  have  a  similar  order  of  magnitude  and  result  in  fairly  uniform 
deposition.  Vertical  surfaces  will  have  a  lower  mass  transfer  rate,  yielding  thicker 
deposits  at  the  edges  of  a  substrate  [28].  Since  the  fibrous  and  porous  substrates  required 
to  make  CMCs  have  both  vertical  and  horizontal  surfaces,  another  method  must  be 
utilized.  At  low  pressures  (-0.001  atm)  the  mass  transfer  rate  becomes  more  ideal  for 
vertical  surfaces.  At  such  low  pressures,  in  order  to  keep  the  same  amount  of  reactants  in 
the  gas  flow,  the  amount  of  carrier  gasses  can  lessened.  Without  the  extra  carrier  gas 
species,  there  is  less  interference  in  adsorption  of  reactants  onto  the  substrate  yielding  in 
higher  deposition  rates  [28],  Increased  temperature  has  a  tendency  to  result  in  smoother, 
more  uniform  deposition  as  well  [29].  This  leads  to  a  common  setup  that  has  an 
inductive  hot  plate  to  heat  the  substrate.  The  walls  of  the  reactor  can  be  heated,  cooled, 
or  insulated.  Hot-wall  reactors  can  aid  in  reducing  the  species  growth  on  surfaces  other 
than  the  substrate  [30]. 
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Silicon  carbide  is  a  highly  variable  material.  It  has  more  than  170  polytypes, 
cubic  (3C  or  (3)  and  hexagonal  (6H  or  a,  and  4H)  being  the  most  common.  It  is  a 
semiconductor,  has  high  hardness,  high  thennal  conductivity,  high  thermal  stability,  and 
high  chemical  resistance  [31].  Further,  as  the  matrix  of  a  CMC  it  becomes  very  fracture 
resistant.  As  silicon  carbide  is  such  a  desirable  ceramic  material,  there  has  been  a  great 
deal  of  research  into  finding  ways  of  producing  it.  Powell  and  Rowland  report  on  a 
laundry  list  of  epitaxy  methods  of  producing  SiC  in  addition  to  CVD  [30].  Even  using 
CVD,  there  are  innumerable  methods  of  producing  SiC.  One  can  use  a  single  precursor 
gas  that  supplies  both  the  silicon  and  the  carbon,  or  one  can  use  two  gasses  that  each 
supply  one  species.  Several  studies  have  been  done  utilizing  methyltrichlorosilane 
(MTS:  CFCSiCb)  as  a  precursor  [6,  10,  12-14,  32-34].  Other  studies  explored  the  use  of 
silicon  tetrachloride  (SiCU)  and  methane  (CH4)  or  other  hydrocarbons  as  a  two-gas 
precursor  system  since  SiCU  is  cheaper  than  MTS  [9,10,  32,  35,  36].  While  the  MTS 
and  SiCU/CFE  systems  were  shown  to  be  effective  in  producing  good  quality  SiC 
deposits,  the  inclusion  of  chlorine  into  the  system  was  found  to  have  a  detrimental  effect 
on  growth  rates.  Chlorine-bearing  species  are  also  harder  to  safely  exhaust  from  the 
reactors  [13,  14].  Removing  the  chlorine  from  the  reaction,  several  studies  have  focused 
on  the  use  of  silane  and  a  simple  hydrocarbon,  usually  methane  or  propane  [5,  7,  8,  10, 

31,  37,  38].  The  relative  low  cost  of  silane  and  propane,  as  well  as  the  lack  of  chlorine  to 
etch  surface  and  cause  problems  when  disposing  of  exhaust  chemicals  led  to  the  choice  of 
using  them  in  this  model  system  to  generate  SiC. 


7 


One  of  the  largest  problems  in  modeling  any  of  these  methods  is  the  sheer  number 
of  reactions  that  occur,  both  in  gas-phase  equilibrium  and  on  the  substrate  surface. 
Stephanie  de  Persis,  et  al.  reported  a  possible  1 1 1  species  with  C-H  bonds  in  610 
irreversible  reactions  and  20  species  with  Si-H  bonds  in  144  irreversible  reactions 
resulting  from  silane  and  hydrogen.  This  did  not  even  include  any  species  with  Si-C 
bonds.  (Chlorinated  species  would  result  in  even  more  species  and  reactions,  yet  another 
good  reason  to  not  use  it  in  an  introductory  model.)  To  reduce  the  computational  power 
needed  for  these  models,  certain  criteria  for  a  species  inclusion  must  be  made.  If  a 
species  is  not  consumed  in  any  reactions  it  can  be  ignored.  If  removing  a  particular 
species  from  the  model  results  in  less  than  a  1%  deviation  in  concentration  profiles,  it  can 
be  ignored  [7].  The  important  reactions  (17  gas-phase  and  15  surface)  and  species  (17) 
chosen  by  Blanquet,  et  al.  are  used  in  this  thesis  as  they  make  up  a  relatively  simple 
model  that  still  has  enough  complexity  to  be  realistic  [5],  Details  of  the  reaction  are  in 
the  section  on  methodology. 

Even  ignoring  redundant  species  and  reactions  to  get  a  reasonably  easy-to- 
compute  model,  all  one  gets  is  information  about  reaction  rates.  To  design  a  true  model 
of  a  deposition  process,  one  must  also  include  the  changes  in  geometry.  William  Ros,  et 
al.  produced  a  method  of  simulating  the  growth  of  a  ceramic  that  relies  on  a  Monte  Carlo 
random  walk  algorithm  [19].  This  is  computationally  expensive  and  focused  only  upon  a 
single  species’  growth.  Ideally,  a  model  should  include  both  detailed  chemical  data  and 
detailed  surface  growth. 
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III.  Modeling  Chemical  Vapor  Deposition 


When  modeling  a  the  chemical  vapor  deposition  process,  there  are  several 
components  that  must  be  included:  bulk  fluid  flow,  diffusion,  chemistry, 
thennodynamics,  and  geometry  changes  resulting  from  fiber  growth.  A  preform  is  placed 
in  a  reactor  and  heated  while  gas-phase  pre-ceramic  chemicals  are  made  to  flow  over  and 
through  the  fibers.  The  chemicals  react  with  each  other,  fonning  new  chemical  species 
that  also  react.  Some  species  will  then  adsorb  to  the  surfaces  of  the  preform.  Once 
adsorbed,  a  species  can  desorb,  stay  on  the  surface  but  remain  inert,  or  react  with  the 
surface  to  create  or  remove  bulk  material.  Those  species  that  desorb  then  continue 
through  the  reactor  and  are  pulled  out  through  the  exhaust  system.  As  surface  reactions 
occur  and  build  up  more  bulk  material,  the  available  paths  of  flow  change  and  can  even 
be  cut  off,  resulting  in  porosity  in  the  bulk  material  [4,  27,  29].  Each  of  these  behaviors 
has  governing  physics  behind  them  that  must  be  understood  in  order  to  model  them 
properly.  This  chapter  discusses  the  background  of  each  area. 

Fluid  Flow 

The  fluid  flow  behavior  of  a  system  is  dependent  on  several  fluid  material 
properties  as  well  as  the  geometry  in  which  the  fluid  moves.  The  motion  of  a  fluid  can  be 
detennined  using  the  thermally-decoupled  Navier-Stokes  Equations.  [Equations  1  &  2] 
They  take  these  particular  forms  since  the  gas-phase  flow  is  compressible.  Equation  1 
represents  the  conservation  of  momentum,  while  Equation  2  shows  the  conservation  of 
mass. 
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As  can  be  seen  in  Equation  1 ,  the  viscosity,  density,  and  pressure  of  a  fluid  are  of 
great  importance  towards  determining  how  it  will  flow  through  a  system.  Increasing 
viscosity  creates  a  more  negative  value  on  the  right  side  of  the  equation,  resulting  in  a 
lower  value  for  the  acceleration  of  the  fluid  (the  time  derivative  of  the  velocity).  A 
similar  results  come  from  increasing  the  density  of  the  fluid..  This  is  to  be  expected,  as 
stickier  and  thicker  fluids  are  harder  to  get  up  to  speed.  The  pressure,  or  more  accurately, 
the  gradient  of  the  pressure  which  can  be  thought  of  as  pressure  difference,  causes  the 
fluid  to  accelerate  towards  the  low  pressure  regions.  And,  since  pressure  is  not  coupled 
to  a  velocity  term,  it  is  sufficient  that  a  pressure  gradient  exists  for  fluid  flow  to  initiate, 
even  without  body  forces  [39,  40], 


Where: 


p^-  +  pu-VU  =  -Vp  +  V-  p(Vu+{Vuf) 


§Av-(ps)=° 


y-(v.„)7  +f 

Equation  1 


Equation  2 


p  =  fluid  density  [kg/m3] 
u  =  velocity  [m/s] 
t  =  time  [s] 
p  =  pressure  [Pa] 
p  =  viscosity  [Pa  s] 

I  =  identity  matrix 
F  =  body  forces  [N/m3] 


The  flowing  that  results  from  pressure  gradients  and  body  forces  can  be,  in 
general,  either  laminar  or  turbulent.  (There  can  be  flows  that  exhibit  both  in  a  mixed 
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region.)  Laminar  flow  is  characterized  by  fluid  moving  in  layers  without  much  cross¬ 
layer  mixing.  Turbulent  flow  is  characterized  by  matter  mixing  between  layers,  rapid 
diffusivity  of  momentum  and  heat,  and  non-zero  vorticity  (curl  of  velocity).  One  method 
that  can  be  useful  in  determining  what  type  of  flow  is  present  in  a  system  is  to  calculate 
the  Reynolds  number  (Re).  It  is  the  ratio  of  inertial  forces  to  viscous  forces.  While  the 
dividing  line  is  often  wide  and  dependent  on  the  geometry,  below  Re  ~  1000  flow  will 
typically  be  laminar;  above  will  typically  be  turbulent.  The  Reynolds  number  is  often 
proportional  to  fluid  velocity,  so  at  the  low  speeds  of  chemical  inflow  that  pre-ceramics 
undergo,  flow  can  be  characterized  as  laminar  [40].  Turbulent  flow  is  usually  considered 
to  be  undesirable  when  using  CVD  as  the  random  flow  tends  to  create  uneven  deposition. 
If  one  can  set  up  a  procedure  where  the  Reynolds  number  remains  in  the  laminar  region, 
deposition  becomes  more  reliable.  According  to  Kern  and  Schnable,  an  ideal  Reynolds 
number  for  atmospheric  pressure  experiments  is  between  10  and  100,  while  for  low 
pressure  setups  it  should  be  closer  to  1  [28], 

Thermodynamics 

Reliable  deposition  is  also  affected  by  the  heat  flowing  through  a  CVD  reactor. 
The  Heat  Equation  [Equation  3],  the  third  Navier-Stokes  equation,  is  used  to  model 
thennal  energy  flow  and  conservation  in  a  system  in  conjunction  with  Equation  2  (mass 
is  still  conserved).  Most  chemical  reaction  rates  in  CVD  are  temperature  dependent  and 
are  endothermic,  changing  the  temperature  of  the  system  further  [29].  The  beauty  of  the 
Heat  Equation  is  that  it  works  in  both  fluids  (the  pre-ceramic  gases)  and  solids  (the 


11 


preform  and  the  reactor  walls).  In  solids,  there  is  no  mass  flow,  so  both  the  velocity  and 
viscous  heating  (Qvh)  terms  are  zero  [39].  In  the  case  of  SiC  CVD,  there  will  not  be 
much  (if  any)  viscous  heating  in  the  gas  phase.  Heat  sources/sinks  include  the  reactor 
wall  temperature,  the  heating  element  for  the  substrate,  and  the  heat  change  that  results 
from  chemical  reaction.  There  is  also  an  element  of  heat  source  and  sink  in  the  inflow 
and  exhaust  of  the  gas  [41].  In  most  CVD  processes,  the  prefonn  substrate  is  heated  by 
either  hot  plate  or  inductive  heating. 


Pcp — + pcPu  ■  vr  =  v  ■  (*vr ) + £ + Qvh 

at 

Equation  3 


Where: 

p  =  density  [kg/m3] 

Cp  =  heat  capacity  [J/kg  K] 

T  =  temperature  [K] 

u  =  fluid  velocity  [m/s] 

k  =  thennal  conductivity  [W/m  K] 

Q  =  heat  sources/sinks  [W/m3] 

Qvh  =  Viscous  heating  source  [W/m3] 


Chemical  Reaction  and  Diffusion 

Chemical  species  move  about  not  only  by  the  convection  applied  by  the  general 
fluid  flow,  but  by  the  diffusion  that  results  from  uneven  concentrations.  This  is  the 
method  by  which  species  adsorb  to  the  surface  of  the  substrate.  The  movement  resulting 
from  diffusion  is  governed  by  Fick’s  Law.  [Equation  4]  The  time  derivative  of  the 
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concentration  accounts  for  accumulation  and  consumption  of  the  species.  The  tenn  with 
the  velocity  component  connects  to  the  convective  transport  from  the  fluid  flow. 

Diffusion  transport  is  represented  by  the  third  tenn.  The  diffusion  coefficient  is  a  tensor 
that  accounts  for  not  only  the  diffusion  of  the  target  species  on  its  own,  but  the  interaction 
it  has  with  each  other  species.  The  final  term  is  used  as  a  source/sink  term,  usually 
related  to  the  combined  rates  of  all  the  chemical  reactions  with  which  the  particular 
species  is  involved  [39]. 


—  +  U-Vc  =  V(DVc)  +  r 

8t 


Equation  4 


Where: 

c  =  concentration  of  a  species  [mol/m3] 
t  =  time  [s] 
u  =  velocity  [m/s] 

D  =  diffusion  coefficient  [m2/s] 
r  =  Reaction  rate  expression  for  species  [mol/  m3  s] 

The  reaction  rate  expression  in  Fick’s  Law  is  a  function  of  the  individual  reactions 
in  which  a  species  occurs.  Specific  reactions  are  governed  by  the  modified  Arrhenius 
equation  for  reaction  rates.  [Equation  5]  It  is  modified  by  the  explicitness  of  the 
temperature  dependence  of  the  pre-factor  A.  Each  of  the  terms  in  the  Arrhenius  equation 
must  be  determined  experimentally  by  measuring  reaction  rates  and  plotting  their  natural 
logarithms  against  the  inverse  of  the  temperature.  The  line  that  results  will  give  -Ea/R  as 
its  slope  and  the  natural  logarithm  of  the  pre  factor  as  the  intercept.  The  temperature 
dependance  factor  P  is  often  viewed  as  a  fudge  factor  to  make  data  fit  a  nice  line,  but  a 
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non-zero  value  can  indicate  that  the  activation  energy  is  ranged.  Individual  species 
diffusion  data  is  detennined  by  the  calculation  of  entropy,  enthalpy,  and  heat  capacity 
from  the  Shomate  Equations  [Equation  6]  [39,  42]. 


r  - 


[x]"[r]m 


Equation  5 


Where: 

r  =  Reaction  rate  [mol/  m3  s] 

[X],  [Y]  =  Concentration  of  reactants  [mol/m3] 
n,  m  =  Reaction  orders 
A=  Pre-factor  [mol 1  -<'n+m)[jn+m)- 1  s' 1  ] 

P  =  Temperature  dependance  factor 
Ea  =  Activation  energy  [J] 

R  =  Gas  constant,  8.3 145  m3  Pa/mol  K 
T  =  Temperature  [K] 

To  =  Reference  temperature,  298  K 


c : 


mol-K 


=  A  +  Bt  +  Ct2  +  Dt 3  +  Et  2 


5° 


J 

_mol-K  _ 


=  At  +  —  Bt2  +  —Ct3  + — Dt4  -Et~l  +F-H 
2  3  4 

=  Alnt  +  Bt  H — Ct~  H — Dt - ~r  +  G 

2  3  2 t2 


Equation  6 


Where: 

Cp°  =  Molar  Heat  Capacity  [J/  mol  K] 

H°  =  Molar  Enthalpy  [kJ/mol] 

S°  =  Molar  Entropy  [J/  mol  K] 
t  =  Temperature  [K]  /  1000K 

A,  B,  C,  D,  E,  F,  G,  H  are  cofactors  determined  by  experiment 
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Deposition  Growth 


Modeling  changes  in  geometry  can  be  difficult.  Random-walk  algorithms 
do  a  good  job,  but  are  computationally  expensive  [19,  26],  A  finite  element  method  can 
be  used  to  model  deposition  growth  by  defonning  the  mesh  itself.  This  can  be 
treacherous,  however.  The  mesh  is  what  the  entire  finite  element  method  is  built  upon.  If 
the  movement  of  the  backbone  has  problems,  the  whole  solution  becomes  problematic. 

The  way  to  use  this  method  is  to  use  a  mesh  with  particularly  small  elements, 
especially  near  the  substrate  surface.  As  mass  adsorbs  onto  the  surface  and  reacts,  the 
gaseous  products  leave  behind  new  surface  species  of  Si  and  C.  Using  knowledge  of  the 
material  density  and  of  the  mass  flow,  retrieved  via  Fick’s  Law,  it  is  possible  to  calculate 
the  velocity  at  which  the  surface  expands  at  each  mesh  element.  This  velocity  can  then 
be  fed  into  the  motion  of  the  mesh  which  indicates  the  growth  of  the  substrate  [39,  43]. 

Finite  Element  Method 

Many  problems  are  very,  very  difficult  to  solve  analytically.  Fortunately,  with  the 
power  of  modern  computing,  it  is  feasible  to  solve  problems  numerically.  The  finite 
element  method  of  numerical  analysis  is  based  upon  the  idea  that  continuous  domains  can 
be  sub-divided  into  finitely  sized,  but  not  necessarily  unifonn,  “elements”.  Over  each 
part  of  the  mesh  of  elements,  an  approximation  of  the  analytical  governing  equations  is 
solved.  These  approximations  are  based  upon  the  integral  forms  of  the  partial  differential 
equations  rather  than  a  Taylor  series  expansion  for  derivatives  as  used  in  the  finite 
difference  method.  Further,  there  is  an  assumption  that  a  solution  can  be  determined  over 
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each  of  these  elements  that  is  a  linear  combination  of  known  functions,  often 
polynomials.  The  individual  solutions  of  the  elements  must  then  be  linked  by  the 
relationships  that  these  elements  have  with  each  other,  usually  by  making  them  into  a 
piece-wise  continuous  solution  over  the  whole  domain.  The  finite  element  method  can, 
however,  handle  discontinuities  with  ease  [44], 

Like  any  numerical  solution,  finite  element  analysis  is  imperfect.  Computers 
must  truncate  the  numbers  stored  in  memory  at  some  physical  limit.  Also,  the  choice  of 
mesh  elements  can  greatly  affect  the  final  result.  If,  for  example,  the  original  domain  to 
be  solved  over  is  a  circle,  it  will  be  impossible  for  a  mesh  of  triangular  elements  to 
completely  cover  the  original  domain.  The  size  of  the  elements  is  also  important.  If  the 
mesh  is  small,  coverage  of  the  original  domain  can  be  made  more  complete  and  the 
solutions  from  the  finite  difference  calculations  will  be  more  accurate.  There  must  be  a 
trade  off  however;  the  smaller  the  mesh  size,  the  more  calculations  must  be  performed 
and  computing  power  has  its  limits  [44].  One  of  the  advantages  that  COMSOL  offers  is 
the  ability  to  customize  mesh  size  in  different  areas  of  the  domain.  Solution  parameters 
often  have  more  complicated  and  abrupt  behaviors  around  edges  and  corners  than  in  bulk 
regions.  A  smaller  mesh  can  be  generated  at  problem  areas  to  get  the  fine  details;  a  larger 
mesh  in  the  bulk  areas  can  help  save  computational  power  [39,  45]. 

IV.  Methodology 

With  an  understanding  of  the  underlying  mechanism  that  govern  the  process  of 
producing  SiC  by  CVD,  the  actual  modeling  process  could  begin.  This  work  was  divided 
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into  two  separate  modeling  systems.  The  first  model  was  intended  as  a  semi-guided 
exploration  into  the  capabilities  of  the  software  used.  A  simple,  supplied  model  was  used 
as  a  starting  point  for  exploring  how  well  the  multiphysics  modules  work  together  and 
what  infonnation  can  be  gleaned  from  such  a  model.  The  second  model  was  aimed 
toward  generating  a  complete,  unified  model  of  both  the  chemistry  and  geometry  changes 
that  occur  when  making  SiC. 

Steam/Propane  Reactor  Tube 

In  an  effort  to  start  understanding  the  physics  involved  with  the  ceramic 
production,  an  example  model  provided  by  COMSOL  was  used  as  a  jumping  point. 

Their  model  for  a  steam  refonner  was  used  for  the  physics  and  chemistry  with  a 
geometry  closer  to  the  CVD  process  to  be  modeled.  This  model  was  of  gaseous  propane 
and  steam  water  reacting  in  a  tube  to  create  carbon  dioxide  and  hydrogen  gas  [46].  A 
complete  report  of  the  model  will  be  included  in  Appendix  D. 

The  model  starts  with  the  two-dimensional  geometry  of  the  reactor  system.  The 
“fibers”  of  the  preform  were  modeled  as  an  array  of  circles  suspended  in  a  space  to  be 
surrounded  by  the  reacting  gasses,  shown  in  Figure  1.  The  offset  arrangement  of  the 
array  of  fibers  was  chosen  to  more  closely  resemble  the  patters  of  woven  fibers  in  a 
fabric.  The  fibers  were  assumed  to  be  incompressible  and  completely  solid  pyrolytic 
carbon  fiber.  The  area  through  which  the  gasses  were  to  flow  was  rectangular, 
representing  a  reactor  sitting  on  a  hot  plate,  having  no  thermal  insulation  on  the  top 
surface,  and  having  an  inlet  and  outlet  that  span  the  complete  vertical  dimension. 


17 


Ambient-FlukJ  Surface 


if 

LL  £ 

Hot  Plate 

Fluid  Flow  Pyrolytic  Carbon  Fibers 

Figure  1,  Steam/propane  reactor:  Geometry  and  setup  of  experiment.  Fiber  radii  are 

0.002  m. 

To  this  geometry,  a  simple  flow  model  was  added  where  material  flows 
from  the  left  side,  the  fluid  inlet,  to  the  right,  the  fluid  outlet.  The  Single  Phase  Laminar 
Flow  module  in  COMSOL  was  appropriate  for  such  flow  as  the  process  of  impregnating 
preforms  is  done  at  low  speeds  and  pressures.  A  no-slip  boundary  condition  was  applied 
to  the  hot  plate,  ambient- fluid  surface,  and  the  fiber  walls.  A  pressure  of  105  Pa  was 
applied  to  both  inlet  and  outlet,  the  inlet  having  an  additional  75  Pa  to  create  a  pressure 
gradient.  The  model  was  run  to  find  the  steady-state  solution.  The  resulting  pressure 
gradient  is  shown  in  Figure  3  with  the  higher  pressure  being  in  red  and  the  lower 
pressures  in  green.  This  pressure  gradient  was  then  the  cause  of  the  fluid  flow  that  is 
shown  in  Figure  4.  As  the  gas  flows  from  high  pressure  to  low  pressure,  the  libers  force 
it  to  flow  into  channels  between  them,  resulting  in  higher  velocity  flow  (shown  in  redder 
tones).  Low  pressure  areas  behind  each  fiber  result  in  small  vortices  of  slower  moving 
(bluer)  fluid. 
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Pressure  Gradient  [PaJ 


Figure  2,  Steam/propane  reactor:  Steady  state  pressure  gradient  [Pa] 
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Figure  3,  Steam/propane  reactor:  Steady  state  velocity  magnitude  [m/s]  and  field  lines. 
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Having  completed  the  flow  portion  of  this  model,  the  thermal  and  chemical 
portions  were  added.  This  utilized  the  Heat  Transfer  in  Fluids,  Heat  Transfer  in  Solids, 
and  Transport  of  Concentrated  Species  modules.  The  starting  temperature  of  the  fluid  in 
the  reactor  and  the  fibers  was  293. 15  K.  The  inflow  fluid  temperature  was  taken  to  be 
700  K  and  the  hot  plate  was  assigned  an  energy  output  of  3000  W/m2.  There  was  heat 
loss  into  the  ambient  atmosphere  on  the  top  surface  and  heat  was  allowed  flow  out  with 
the  fluid.  The  reaction  C3H8  +  H2O  ->  H2  +  CCkwas  entered  into  the  system  along  with 
initial  mass  fractions  for  each  species  (0.28,  0.50,  0.01,  and  0.01  respectively).  The  heat 
of  the  reaction  was  also  included  into  the  Heat  Transfer  module.  These  updates  to  the 
model  resulted  in  the  steady-state  solutions  for  temperature  and  the  mass-fraction  of 
propane  shown  in  Figures  4  and  5.  Figure  4  shows  a  temperature  gradient  with  red  being 
the  hottest  and  blue  being  the  coolest.  As  could  be  expected,  the  cooler  temperatures 
were  along  the  ambient-fluid  boundary  where  heat  is  lost  as  well  as  towards  the  outlet. 
This  also  makes  sense  as  the  reaction  was  endothennic  and  had  a  fairly  high  activation 
energy  resulting  in  much  of  the  heat  energy  available  being  used.  This  also  corresponded 
to  the  mass  fraction  of  propane  map  in  Figure  5.  At  the  onset  (near  the  inlet)  there  was  a 
higher  ratio  of  propane  in  the  system  as  its  had  not  has  the  opportunity  of  react.  There 
was  also  a  higher  ratio  at  the  fluid-ambient  surface  as  there  was  not  as  much  of  the  heat 
that  was  needed  to  activate  the  reaction.  The  low  pressure  vortices  seen  in  Figure  3  held 
the  fluid  for  more  time,  giving  the  chemicals  more  chance  to  react.  This  lead  to  the  lower 
mass  fraction  of  propane  behind  each  fiber.  Similar  data  for  the  other  chemical  species, 
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Figure  4,  Steam/propane  reactor:  Steady  state  temperature  [K] 
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Figure  5,  Steam/propane  reactor:  Steady  state  mass  fraction  of  C3H8 
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though  not  shown,  can  be  drawn  out  the  model  with  ease.  In  a  SiC  deposition  system, 
these  type  of  features  in  the  flow  may  cause  nonuniform  deposition  on  the  substrates. 
This  shows  how  important  it  is  to  optimize  the  system  parameters. 

This  initial  exploration  into  the  capability  that  COMSOL  offered  was  quite 
revealing.  COMSOL  was  able  to  track  a  great  many  physical  characteristics  (i.e. 
temperature,  velocity,  and  concentration)  simultaneously  and  was  able  to  report  back  on 
each  at  any  given  point  or  average  across  a  surface.  This  shows  good  promise  for  being 
able  to  generate  the  desired  complete  SiC  model. 

Silicon  Carbide  Model  System  in  2-D 

After  exploring  the  simple  propane  and  steam  system,  the  more  complex 
silicon  carbide  deposition  process  was  tackled.  The  new  phase  of  exploration  was  done 
in  three  disparate  parts:  a  flow  model,  a  chemical  reaction  model,  and  a  moving  mesh 
model.  With  each  part,  an  overly-simplified  model  was  used  in  order  to  understand  the 
way  that  COMSOL  models  the  physics  and  what  kinds  of  effects  various  changes  in  the 
initial  settings  can  have. 

Flow  Model 

The  flow  model  started  out  with  a  very  similar  geometry  to  the  initial 
exploration  model.  An  off-set  array  of  circular  “fibers”  around  which  gas-phase  fluid 
flows  was  generated  as  in  Figure  6.  Each  fiber  was  taken  to  have  a  radius  of  0.002  m.  A 
rectangular  area  surrounds  the  fibers  as  a  model  of  the  reactor.  An  additional  rectangle 
was  placed  in  between  the  inlet  and  the  fibers  to  start  including  larger,  more  complex 
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geometries.  The  Single-Phase  Laminar  Flow  module  was  used.  The  inlet  is  the  left-most 
edge  and  the  outlet,  the  right-most.  A  fluid  pressure  difference  was  applied  between  inlet 
and  outlet  to  cause  flow.  At  the  onset,  the  same  fluid  properties  were  used  as  in  the 
previous  model.  Thermodynamics  were  ignored  for  this  portion  of  the  model,  so  there 
exists  no  hot  plate  or  ambient  heat  loss. 

Velocity  Magnitude  [m  s) 
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Figure  6,  Flow  model:  Geometry  (distances  in  [m])  and  velocity  magnitude  (color 
gradient)  [m/s].  Pressure  difference  of  491  Pa  and  X,  Y  fiber  distance  of  thrice  the  radius. 

Since  the  introductory  exploration  of  steam  and  propane  was  only  looking  at  a 
static  model,  and  understanding  how  optimization  of  said  model  works  is  part  of  the  goal 
in  this  endeavor,  the  parameterization  of  variables  in  the  model  was  investigate.  In  this 
flow  portion  of  the  model,  there  are  only  a  few  things  that  can  be  parameterized:  the 
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pressure  difference  between  inlet  and  outlet,  the  geometry,  and  the  fluid  properties.  The 
fluid  properties  are  to  be  controlled  by  the  chemistry  and  thermodynamics,  so  the 
pressure  and  geometry  were  focused  on.  In  this  first  case,  given  a  geometry  with  each 
fiber  being  separated  in  X  and  Y  direction  by  thrice  the  fiber  radius,  the  pressure 
difference  was  parameterized  from  0. 1  to  50  Pa.  The  average  fluid  velocity  taken  across 
the  outlet  edge  was  calculated  [Figure  7].  The  maximum  Reynolds  number  taken  from 
every  finite  element  mesh  cell  was  also  calculated  [Figure  8],  As  long  as  the  fluid 
remains  in  the  laminar  flow  region,  the  deposition  should  be  fairly  stable.  The  Reynolds 
number  is  useful  in  helping  to  detennine  how  turbulent  the  flow  was. 


Average  Outlet  Velocity  v.  Pressure  Difference 


Figure  7,  Flow  model:  Average  outlet  velocity  [m/s]  for  pressure  differences  from  0.1  to 
50  Pa  and  X,Y  separation  of  thrice  fiber  radius 
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As  the  geometries  used  in  this  thesis  were  deliberately  simplified,  it  was  not  of 
great  importance  to  investigate  the  effects  of  differences  of  geometry  in  detail,  but  in  the 
future,  such  investigation  could  be  most  helpful.  One  example  was  investigated: 
changing  the  X-directional  separation  of  the  fibers  with  a  constant  difference  of  inflow 
pressure  to  outflow  pressure  of  3  Pa.  The  separation  ranged  from  twice  the  fiber  radius 
(fibers  almost  touching)  to  twenty  times  the  radius.  This  modeled  the  differences 
between  tighter  and  looser  weaves  of  fabric.  As  with  pressure  differences,  average  outlet 
velocity  and  maximum  Reynolds  number  were  calculated  [Figures  9  and  10]. 

Maximum  Cdl  Reynolds  Number  v.  Pressure  Difference 


Figure  8,  Flow  model:  Maximum  Reynolds  number  for  pressure  differences  from  0.1  to 
50  Pa  and  X,  Y  separation  of  thrice  fiber  radius. 

Figure  9  shows  that  a  maximum  average  output  flow  rate  can  be  found  when  the 

horizontal  separation  of  fibers  is  between  0.006  and  0.008  m.  Since  in  this  model  the 

fibers  are  of  radius  0.002  m,  the  maximum  average  output  flow  rate  is  found  when  the 
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fiber  centers  are  separated  at  a  distance  between  3  and  4  times  the  fiber  radius.  While  the 
Reynolds  number  seems  to  be  somewhat  random  at  first  look,  it  only  ranges  from  9.4  to 
12.2.  When  the  pressure  was  varied,  the  Reynolds  number  ranged  from  0  to  more  than  40 
[Figure  8].  This  would  suggest  that  the  tightness  of  fabric  weave  does  not  have  as  great 
an  effect  on  fluid  turbulence  as  the  pressure  difference  from  inlet  to  outlet  does.  This 
would  further  indicate  that  while  the  choice  of  fabric  style  will  matter  to  the  final  material 


characteristics,  it  will  not  have  great  effect  on  the  production  characteristics. 


Average  Outlet  Velocity  v.  x-Uistanec  Between  fibers 


Average 

Outlet 
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x-L>istance  Between  Centers  of  I ibers  [f  iber  Radius | 

Figure  9,  Flow  model:  Average  outlet  velocity  for  fiber  x-distance  from  2-20  fiber  radii 

and  pressure  difference  of  3  Pa. 
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x-Distancc  Between  Centers  of  Fibers  (Fiber  Radius | 

Figure  10,  Flow  model:  Maximum  Reynolds  number  for  fiber  x-distance  from  2-20  fiber 

radii  and  pressure  difference  of  3  Pa. 

Chemical  Reaction  Model 

The  second  model  required  for  this  project  was  the  of  chemical  reactions 
and  their  associated  thermodynamics.  The  process  of  converting  silane  and  propane  to 
silicon  carbide  in  a  hydrogen  environment  is  complex  and  involves  well  over  600 
reactions  [7,  8],  A  reduced  model  was  used  as  reported  by  Blanquet,  et  al.  with  32 
reactions  shown  in  Tables  1  and  2[5],  The  first  two  reactions  in  Table  2  indicate  that 
hydrogen  can  etch  the  surface.  Reaction  1 1  in  Table  1  is  incompletely  reported.  Since 
Reaction  12  from  Table  1  and  Reaction  10  from  Table  2  are  both  tied  to  it,  and  both 
silicon  and  carbon  are  deposited  without  them,  all  three  were  ignored  in  this  model 
building  process.  This  was  not  thought  to  be  too  egregious  as  Dollet  reported  that 
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H3S1CH3  and  similar  species  were  not  created  in  large  amounts  (at  1573  K  and  5332  Pa) 
[8].  It  was  reported  that  SbC  and  C2H2  were  the  main  gas  species  that  were  involved 
with  the  deposition  of  SiC  directly  when  performed  at  a  pressure  of  250  mbar  and 
temperatures  from  1600  to  2000°C.  At  higher  temperature  and  lower  pressure,  SiC2 
replaces  C2H2  as  one  of  the  top  contributing  elements  [5]. 

The  Chemical  Reaction  Engineering  module  was  used  in  a  non-dimensional, 
time-dependent  study.  The  fifteen  reactions  of  the  gas  phase  were  entered  into  the  model 
in  the  form  of  coefficients  of  the  rate  equation  [Equation  5].  The  coefficients  of  the 
Shomate  Equations  [Equation  6]  for  each  species  were  also  entered.  While  Blanquet,  et 
al.  reported  the  forward  reaction  rate  coefficients  they  used  in  their  models,  there  were  no 
reverse  reaction  coefficients  [5].  It  was  decided  to  search  for  all  coefficients  to  all 
equations  in  the  Chemistry  WebBook  posted  by  the  National  Institute  for  Standards  and 
Technology  (NIST)  [42],  These  values  are  reported  in  Appendix  A  for  the  species  data 
and  Appendix  B  for  the  reaction  data.  The  initial  ratio  of  Silane/Propane/Hydrogen  was 
3/1/8.  Running  the  model  for  from  (0  to  10)*  10'7  s  at  a  temperature  of  2073  K  resulted  in 
the  concentrations  of  gaseous  species  shown  in  Figure  11.  Figure  12  shows  the  same,  but 
zoomed  in  to  see  the  lower  concentration  species.  Hydrogen  gas,  H2  remains  near  a 
concentration  of  0.5  mol/m3,  so  it  was  not  shown  on  the  figures.  Silane  comes  to  a  steady 
state  of  about  0. 12  mol/m3,  while  the  propane  is  used  up  entirely.  Hydrogen  atoms, 
methane,  silylene,  and  atomic  silicon  make  up  the  bulk  of  the  remaining  gas  phase 
species. 
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Table  1,  Gas  phase  reactions  and  their  Arrhenius  equation  parameters  [5], 
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Table  2,  Surface  reactions  and  their  Arrhenius  equation  parameters  [5]. 


30 


There  was  difficulty  in  getting  the  Reaction  Engineering  module  to  work 
smoothly.  The  built-in  versions  for  the  Shomate  equations  [Equation  6]  were  different 
than  the  ones  reported  by  NIST.  The  cofactors,  thus,  could  not  be  entered  in  without 
errors.  A  work-around  was  attempted,  but  time  constraints  prevented  the  finished  model 
from  being  built.  Had  this  been  completed,  experiments  would  have  been  done  where 


temperature  would  be  parameterized  as  well  as  initial  concentrations  of  species. 


Girscous  Species  Conecntralion  v.  Time 


Figure  11:  Chemistry  model:  Molar  concentrations  [mol/m3]  of  gas-phase  species  as  a 

function  of  time  [10'7  s]. 
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Gaseous  Species  Concentration  \.  Time 


Figure  12,  Chemistry  model:  Molar  concentrations  [mol/m3]  of  gas-phase  species  as  a 

function  of  time  [10'7  s],  zoomed  in. 


Moving  Mesh 

The  next  modeling  step  was  to  create  a  model  that  would  include  the  deformation 
of  the  geometry  of  the  system  as  chemical  reactions  deposit  more  silicon  carbide  on  the 
surfaces.  A  combination  of  the  Transport  of  Diluted  Species  and  Deformed  Geometry 
modules  was  used.  To  give  a  proof  of  concept,  rather  than  the  complex  chemistry  from 
the  previous  model,  as  simple  single  species  model  was  created.  The  geometry  was  of  a 
square  of  0.1  m  for  the  matrix  area  and  four  circular  fibers  of  radius  0.005  m,  shown  in 
Figure  13.  An  ideal  gas  with  an  initial  concentration  of  20  mol/m3,  a  diffusion  coefficient 
of  10'9  m2/s,  and  a  density  of  10  mol/m3  when  in  solid  form  was  used.  The  concentration 
of  the  gas  phase  was  contained  to  be  always  zero  at  the  boundaries  of  the  fibers.  In  a 
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more  true-to-life  model,  the  concentration  of  the  species  at  the  boundaries  would  be 
determined  by  the  chemical  reaction  rates  that  produce  SiC  using  Fick’s  Law. 

The  mesh  used  for  the  finite  element  analysis  has  to  have  constraints  put  upon  it 
for  the  deformation  to  occur.  The  mesh  of  the  matrix  area  was  given  the  freedom  to 
move  without  any  constraint.  The  mesh  along  outer  boundaries  on  left  and  right  edges 
was  constrained  to  only  move  in  vertical  directions.  Similarly,  the  top  and  bottom  edges 
had  a  mesh  constrained  to  only  move  horizontally.  Edges,  after  all,  cannot  move  outward 
in  an  enclosed  volume.  The  boundaries  of  the  libers  were  constrained  to  a  velocity  that 
results  from  flow  rates  as  detennined  by  Fick’s  Law  for  diffusion  and  by  material  density. 


0.1 

0.08 

0.06 
0.04 

0.02 
0 

Figure  13,  Moving  mesh  model:  Geometry.  Sides  are  0.1  m,  circle  radii  are  0.005  m. 
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The  model  was  then  run  as  a  time  dependent  study  in  steps  of  100  s  for  106  s.  The 
results  are  shown  in  Figure  14.  The  white  circles  are  the  original  fiber  positions  and 
sizes.  The  grey  areas  being  the  accumulated  bulk  material.  The  colored  areas  show  the 
concentration  of  the  remaining  gas  species.  One  problem  with  this  method  is  that  the 
mesh  tends  to  become  so  defonned  that  the  mesh  elements  will  invert  and  start  modeling 
unreal  physics.  To  combat  this,  the  model  has  to  occasionally  be  stopped,  a  new  mesh  set 
to  the  new  geometry,  and  started  up  again.  An  automatic  re-meshing  can  be  set  up,  but  it 
tends  to  keep  odd  artifacts  from  the  previous  mesh. 

Concentration  of  Gas  Species  [mul/tnA3] 

0.1 


0.05 


0 

m  0  0.05  0.1 

Figure  14,  Moving  mesh  model:  Surface  growth  model  after  time  lapse.  White  is  original 
geometry,  grey  deposited  material,  color  scale  remaining  species  concentration  [mol/m3]. 
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V.  Conclusions  and  Recommendations 


Conclusions  of  Research 

This  research  project  investigated  the  creation  of  a  series  of  models  that  would 
build  towards  a  unified  chemical  and  geometric  model  of  a  SiC  CVD  process.  This  was 
done  to  gain  a  better  understanding  of  how  the  CVD  process  works  and  how  the  given 
parameters  affect  the  final  product.  This  was  also  done  in  an  effort  to  find  a  modeling 
software  package  that  can  handle  rigors  of  modeling  the  complex  process  that  results  in 
such  desirable  materials. 

The  first  model,  that  of  a  steam  and  propane  reactor  tube,  showed  how  well 
COMSOL  software  is  at  pulling  out  material  information  from  the  model  at  all  times  and 
positions.  Velocity,  pressure,  mass  fraction,  and  temperature  were  all  very  easy  to  attain. 
The  second  model  helped  to  showed  some  of  the  effects  of  changes  in  pressure  and 
geometry.  The  third  model  was  unfortunately  not  completed  due  to  time  constraints.  It 
was,  however,  able  to  at  least  verify  that  propane,  as  the  carbon  source  and  limiting  factor 
in  deposition  rate,  is  consumed  very  quickly  [33].  The  final  model  showed  the  feasibility 
of  using  a  finite  element  method  for  modeling  deposition  growth  without  the  expense  that 
comes  from  random  walk  algorithms. 

Significance  of  Research 

This  thesis  was  a  preliminary  study  into  using  COMSOL  as  a  method  for 
modelingSiC  CVD  production.  It  showed  that  the  software  is  capable,  even  if 
occasionally  problematic  (like  with  the  Shomate  equation  parameters).  The  various 
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parameters  that  are  given  as  input  into  the  model  have  different  levels  of  effect  on  the 
resulting  simulations.  Fluid  pressure  is  a  major  component  of  determining  how  uniform 
the  flow  is.  Lower  pressure  results  in  more  uniform  (laminar)  flow.  Even  with  laminar 
flow,  vortices  of  low  pressure  can  form.  This  leads  to  uneven  reaction  rates  and  species 
concentrations  across  the  substrate  surfaces.  This  can  strongly  affect  the  uniformity  of 
the  deposition  of  SiC  on  the  substrate.  Temperature  has  a  large  effect  on  deposition  rates, 
but  as  the  chemical  reaction  model  was  unable  to  be  finished,  it  was  unable  to  be  verified 
numerically.  Changes  in  initial  geometry  did  not  have  much  effect  on  flow 
characteristics. 

Recommendations  for  Future  Research 

A  completely  integrated  model  is  within  reach,  it  was  only  time  constraints  that 
prevent  its  inclusion  in  this  thesis.  Once  a  complete  model  for  SiC/substrate  production 
via  CVD  is  complete,  it  would  be  advantageous  to  compare  the  model  to  actual 
experiment  to  see  how  well  it  correlates  with  reality  and  what  changes  might  need  to  be 
made.  Additionally,  more  complex  (and  therefor  realistic)  geometries  should  be  used. 
(An  example  three-dimensional  geometry  is  included  in  Appendix  C.) 

It  is  possible  for  users  to  build  their  own  modules  for  use  in  COMSOL.  If  the 
time  and  expertise  can  be  applied  to  the  work,  a  specialized  module  for  modeling  the 
CVD  process  could  be  created.  Furthermore,  it  would  be  entirely  feasible  to  generate 
similar  modules  for  other  forms  of  CMC  production,  such  as  polymer  liquid  infiltration. 
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Appendix  A:  Chemical  Species  Data 


This  section  lists  the  data  of  individual  species  as  taken  from  the  NIST  Chemistry 
WebBook.  Species  name,  formula,  and  molecular  mass  is  given  on  the  first  line  for  each 
species.  The  ranges  of  temperatures  (in  Kelvin)  over  which  the  Shomate  parameters  (A- 
H)  are  valid  are  given  in  the  second  line.  Shomate  parameter  units  are  as  follows:  A  [J/ 
mol  K],  B  [J/mol  K2],  C  [J/mol  K3],  D  [J/mol  K4],  E  [J  K/mol],  F  [J/mol],  G  [J/mol  K],  H 
[J/mol],  [30] 

Table  3,  Properties  of  individual  chemical  species.  Includes  name,  formula,  molecular 
mass,  and  Shomate  parameters  with  the  temperatures  over  which  they  are  valid. 


Silane 

SiH 

32.1173 

298.-  1300.  K 

1300.  -  6000.  K 

A  [J/mol  K] 

6.060189 

99.84949 

B  [J/mol  K 

139.9632 

4.251530 

C  [J/mol  K 

-77.88474 

-0.809269 

D  [J/mol  K 

16.24095 

0.053437 

E  [J  K/mol] 

0.135509 

-20.39005 

F  [J/mol] 

27.39081 

-40.54016 

G  [J/mol  K] 

174.3351 

266.8015 

H  [J/mol] 

34.30905 

34.30905 

Silylene 

:SiH2 

30.1014 

A  [J/mol  K] 

B  [J/mol  K 

C  [J/mol  K 

D  [J/mol  K 

E  [J  K/mol] 
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F  [J/mol] 

G  [J/mol  K] 

H  [J/mol] 

Hydrogen 

H 

2.01588 

Gas 

298.  -  1000  K. 

1000. -2500.  K 

2500.  -  6000.  K 

A  [J/mol  K] 

33.066178 

18.563083 

43.413560 

B  [J/mol  K 

-11.363417 

12.257357 

-4.293079 

C  [J/mol  K 

11.432816 

-2.859786 

1.272428 

D  [J/mol  K 

-2.772874 

0.268238 

-0.096876 

E  [J  K/mol] 

-0.158558 

1.977990 

-20.533862 

F  [J/mol] 

-9.980797 

-1.147438 

-38.515158 

G  [J/mol  K] 

172.707974 

156.288133 

162.081354 

H  [J/mol] 

0.0 

0.0 

0.0 

Disilane 

Si 

62.2186 

A 

[J/mol  K] 

B 

[J/mol  K 

C 

[J/mol  K 

D 

[J/mol  K 

E 

[J  K/mol] 

F 

[J/mol] 

G 

[J/mol  K] 

H 

[J/mol] 

Silicon 

Si 

28.0855 

298.  -  1685  K. 

1685.-  3505.  K 

3505.-6000.  K 

A 

[J/mol  K] 

22.81719 

27.19604 

14.59321 

B 

[J/mol  K 

3.899510 

-1.198306x10 

5.224644 

C 

[J/mol  K 

-0.082885 

5.353262x10 

-1.078879 

D 

[J/mol  K 

0.042111 

-6.956612x10 

0.074000 

E 

[J  K/mol] 

-0.354063 

-4.294375x10 

2.309405 

450.3365 

190.2494 

450.0018 


F  [J/mol] 

-8.163946 

40.36163 

450.3365 

G  [J/mol  K] 

43.27846 

77.37178 

190.2494 

H  [J/mol] 

0.000000 

48.46997 

450.0018 

Hydrogen 

*H 

1.00794 

Elemental 

298.  -  6000.  K 

A  [J/mol  K] 

20.78603 

B  [J/mol  K 

4.850638x10 

C  [J/mol  K 

-1.582916x10 

D  [J/mol  K 

1.525102x10 

E  [J  K/mol] 

3.196347x10 

F  [J/mol] 

211.8020 

G  [J/mol  K] 

139.8711 

H  [J/mol] 

217.9994 

Propane 

C 

44.0956 

A  [J/mol  K] 

B  [J/mol  K 

C  [J/mol  K 

D  [J/mol  K 

E  [J  K/mol] 

F  [J/mol] 

G  [J/mol  K] 

H  [J/mol] 

Methyl  Radical 

*CH 

15.0345 

298.  -  1400.  K 

1400.  -  6000.  K 

A  [J/mol  K] 

28.13786 

67.18081 

B  [J/mol  K 

36.74736 

7.846423 

C  [J/mol  K 

-4.347218 

-1.440899 

D  [J/mol  K 

-1.595673 

0.092685 

E  [J  K/mol] 

0.001860 

-17.66133 

A  [J/mol  K] 
B  [J/mol  K 
C  [J/mol  K 
D  [J/mol  K 
E  [J  K/mol] 
F  [J/mol] 

G  [J/mol  K] 
H  [J/mol] 


A  [J/mol  K] 
B  [J/mol  K 
C  [J/mol  K 
D  [J/mol  K 
E  [J  K/mol] 
F  [J/mol] 

G  [J/mol  K] 
H  [J/mol] 


A  [J/mol  K] 
B  [J/mol  K 
C  [J/mol  K 
D  [J/mol  K 
E  [J  K/mol] 
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F  [J/mol] 

G  [J/mol  K] 

H  [J/mol] 

Ethylene 

C 

28.0532 

298.  -  1200.  K 

1200.  -  6000.  K 

A  [J/mol  K] 

-6.387880 

106.5104 

B  [J/mol  K 

184.4019 

13.73260 

C  [J/mol  K 

-112.9718 

-2.628481 

D  [J/mol  K 

28.49593 

0.174595 

E  [J  K/mol] 

0.315540 

-26.14469 

F  [J/mol] 

48.17332 

-35.36237 

G  [J/mol  K] 

163.1568 

275.0424 

H  [J/mol] 

52.46694 

52.46694 

Acetylene 

C 

26.0373 

A  [J/mol  K] 

B  [J/mol  K 

C  [J/mol  K 

D  [J/mol  K 

E  [J  K/mol] 

F  [J/mol] 

G  [J/mol  K] 

H  [J/mol] 

Silicon  dimer 

*Si 

56.171 

298.-  1100.  K 

1100. -6000.  K 

A  [J/mol  K] 

24.19988 

47.92605 

B  [J/mol  K 

31.58280 

-4.345921 

C  [J/mol  K 

-8.262313 

1.037265 

D  [J/mol  K 

-3.875582 

-0.072474 

E  [J  K/mol] 

0.150614 

-0.208383 

F  [J/mol] 

581.9149 

571.7143 

G  [J/mol  K] 

250.9182 

280.7259 

H  [J/mol] 

589.9440 

589.9440 

Disilicon  Carbide 

Si 

68.1817 

298.  -  1200.  K 

1200.  -  6000.  K 

A  [J/mol  K] 

38.96220 

62.75874 

B  [J/mol  K 

41.84126 

-7.005439 

C  [J/mol  K 

-29.86230 

4.097706 

D  [J/mol  K 

7.853285 

-0.425546 

E  [J  K/mol] 

-0.139347 

0.696590 

F  [J/mol] 

521.8578 

513.9207 

G  [J/mol  K] 

277.4360 

312.3410 

H  [J/mol] 

535.5520 

535.5520 

SiCH 

42.1121 

A  [J/mol  K] 

B  [J/mol  K 

C  [J/mol  K 

D  [J/mol  K 

E  [J  K/mol] 

F  [J/mol] 

G  [J/mol  K] 

H  [J/mol] 

Methylene 

:CH 

14.0266 

298.  -  1400.  K 

1400.  -  6000.  K 

A  [J/mol  K] 

31.96823 

51.55901 

B  [J/mol  K 

6.783603 

3.876975 

C  [J/mol  K 

12.51890 

-0.649608 

D  [J/mol  K 

-5.696265 

0.037901 

E  [J  K/mol] 

-0.031115 

-10.72589 

F  [J/mol] 

376.3558 

350.6715 

G  [J/mol  K] 

229.9150 

232.3212 

H  [J/mol] 

386.3924 

386.3924 

Carbon 

C 

12.0107 

Elemental 

298.  -  6000.  K 

A  [J/mol  K] 

21.17510 

B  [J/mol  K 

-0.812428 

C  [J/mol  K 

0.448537 

D  [J/mol  K 

-0.043256 

E  [J  K/mol] 

-0.013103 

F  [J/mol] 

710.3470 

G  [J/mol  K] 

183.8734 

H  [J/mol] 

716.6690 

Silicon  Carbide 

SiC 

40.0962 

Gas  Phase 

298.  -  1000.  K 

1000.  -  6000.  K 

A  [J/mol  K] 

60.32180 

33.36570 

B  [J/mol  K 

4.423940 

4.437140 

C  [J/mol  K 

-61.79470 

-0.865113 

D  [J/mol  K 

39.06140 

0.057726 

E  [J  K/mol] 

-1.675991 

3.073291 

F  [J/mol] 

696.3070 

717.4580 

G  [J/mol  K] 

277.6700 

262.5750 

H  [J/mol] 

719.6480 

719.6480 

Silicon  Carbide  alpha 

298.  -  1000.  K 

1000. -4000.  K 

A  [J/mol  K] 

20.55859 

46.90222 

B  [J/mol  K 

64.57962 

5.845968 

C  [J/mol  K 

-52.98827 

-1.085410 

D  [J/mol  K 

16.95813 

0.093021 

E  [J  K/mol] 

-0.781847 

-3.448876 

F  [J/mol] 

-82.73693 

-95.46716 
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G  [J/mol  K] 

19.90848 

56.97520 

H  [J/mol] 

-71.54598 

-71.54598 

Silicon  Carbide  beta 

298.  -  1100. 

1100.-4000. 

A  [J/mol  K] 

20.50009 

48.22227 

B  [J/mol  K 

63.37170 

5.004148 

C  [J/mol  K 

-49.54023 

-1.037594 

D  [J/mol  K 

14.82801 

0.086339 

E  [J  K/mol] 

-0.759969 

-3.912333 

F  [J/mol] 

-84.29337 

-98.46542 

G  [J/mol  K] 

20.30926 

57.76221 

H  [J/mol] 

-73.22000 

-73.22000 
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Appendix  B:  Chemical  Reaction  Data 


The  following  tables  list  the  coefficients  for  the  Arrhenius  equations  (Equation  #) 
for  each  homogenous  reaction  supplied  by  Blanquet,  et  al.  (Table  1)  The  coefficients 
were  taken  from  the  NIST  Chemistry  WebBook.  The  units  for  the  reaction  rate 
coefficient  k  change  with  the  order  of  the  reaction,  thus  the  units  for  the  Arrhenius  pre¬ 
factor  (A)  change  as  well.  A  is  in  units  of  [s1]  for  reactions  of  order  1,  [cm3  mol'1  s_1]  for 
reactions  of  order  2,  and  [cm6  mol'2  s'1]  for  reactions  of  order  3.  For  those  reactions 
where  NIST  was  not  able  to  supply  data,  the  forward  reaction  cofactors  were  used  from 
Blanquet,  et  al.  The  reverse  reaction  cofactors  were  calculated  based  upon  the 
thermochemistry  of  the  individual  species.  [2,  30] 


Table  4,  Homogeneous  reaction  forward  reaction  rate  equation  parameters.  Reaction 
numbers  correspond  to  those  in  Table  1 .  Values  retrieved  from  NIST  Chemical  Kinetics 
Database.  Units  of  A  are  [s'1]  for  reactions  of  order  1,  [cm3  mol'1  s'1]  for  reactions  of 
order  2,  and  [cm6  mol'2  s'1]  for  reactions  of  order  3.  The  reference  temperature  is  298  K. 


Forward 

Reaction 

Reaction  Order 

A 

n 

Ea  [J] 

1 

1 

4.92E+28 

-15.9 

3.42E+04 

2 

1.67E-06 

18.2 

-6.26E+04 

2 

1 

8.96E+22 

-9.7 

2.62E+05 

3 

1 

9.73E+08 

0.0 

2.61E+05 

2 

3.90E+28 

-12.5 

2.87E+05 

4 

5 

1 

5.42E+18 

-2.3 

3.71E+05 
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2 

2.26E+30 

-11.2 

3.66E+05 

6 

2 

6.56E+11 

2.7 

3.73E+04 

7 

2 

2.82E+13 

0.1 

-3.27E+02 

8 

2 

3.27E+13 

-0.5 

-1.67E+02 

3 

1.40E+24 

-7.1 

1.32E+04 

9 

2 

2.43E+13 

0.1 

9.03E+03 

3 

5.64E+18 

-0.2 

4.65E+03 

10 

1 

7.50E+06 

5.0 

2.41E+05 

2 

3.61E+24 

-6.2 

4.08E+05 

11 

12 

13 

14 

15 

16 

17 

Table  5,  Homogeneous  reaction  reverse  reaction  rate  equation  parameters.  Reaction 
numbers  correspond  to  those  in  Table  1 .  Values  retrieved  from  NIST  Chemical  Kinetics 
Database.  Units  of  A  are  [s'1]  for  reactions  of  order  1,  [cm3  mol'1  s'1]  for  reactions  of 
order  2,  and  [cm6  mol'2  s'1]  for  reactions  of  order  3.  The  reference  temperature  is  298  K. 


Reverse 

Reaction 

Reaction  Order 

A 

n 

Ea  [J] 

1 

2 

2.17E+08 

3.1 

-1.06E+04 

2 

2 

3.13E+14 

-1.1 

3.45E+02 

3 

4 
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5 

2 

1.72E+13 

-0.2 

-9.81E+02 

3 

4.25E+27 

-16.1 

1.58E+04 

6 

2 

2.46E+11 

1.6 

4.22E+04 

7 

2 

4.03E+12 

0.6 

4.26E+04 

8 

1 

6.92E+18 

-2.7 

3.86E+05 

2 

4.28E+17 

1.5 

3.78E+05 

9 

1 

5.85E+12 

0.6 

1.60E+05 

2 

7.06E+18 

-0.8 

1.54E+05 

10 

2 

3.01E+11 

0.0 

1.63E+05 

11 

12 

13 

3 

5.94E+14 

0.0 

1.97E+03 

14 

15 

16 

17 
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Appendix  C:  3-D  Fabric  Preform  Model 


After  completing  the  two  dimensional  model,  it  would  be  considered  ideal 
to  extend  into  three  dimensions.  This  would  result  in  a  more  realistic  geometry,  and  thus 
a  more  realistic  model.  The  geometry  of  this  models  was  chosen  to  be  layers  of  simple 
square- woven  fabric. 

Modeling  the  fabric  started  with  a  single  fiber.  It  was  modeled  as  the 
sweep  of  a  circular  face  along  a  sinusoidal  backbone.  The  fiber  was  then  replicated 
thrice.  Each  replica  was  then  translated,  flipped,  and/or  rotated  until  there  were  four 
fiber,  two  in  the  warp  and  two  in  the  weft,  in  the  beginning  of  a  square -weave.  The  warp 
and  weft  fibers  were  then  replicated  horizontally  until  the  desired  size  of  the  fabric  was 
reached.  Then,  the  resulting  single  layer  of  fabric  was  then  replicated  vertically  to  create 
the  complete  preform  as  shown  in  Figure  15.  A  COMSOL  model  report  follows. 


Figure  15,  3-D  fabric  preform  model. 
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3-D  Fabric  Geometry 


1.  Global  Definitions 
1.1.  Parameters  1 
Parameters 


Name 

Expression 

Description 

R  fib 

2 

Radius  of  fiber 

nx_fib 

10 

number  of  fibers  in  x  direction 

n  layer 

5 

number  of  layers 

ny_fib 

12 

number  of  fibers  in  y  direction 
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3-D  Fabric  Geometry 


2.  Model  1  (modi) 

2.1.  Definitions 

2.1.1.  Selections 

Fibers 

Selection  type 

Explicit 

I - 

Selection 

Domains  2-111 


Fiber  Boundaries 

Selection  type 

Explicit 

Selection 

Boundaries  6-665 


2.1.2.  Coordinate  Systems 
Boundary  System  1 


Coordinate  system  type 

Boundary  system 

Identifier 

sysl 

Settings 


Name 

Value 

Coordinate  names 

{tl,t2,  n} 

Create  first  tangent  direction  from 

Global  Cartesian 
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3-D  Fabric  Geometry 


2.2.  Geometry  1 


Geometry  1 

Units 


Length  unit 

m 

Angular  unit 

deg 

Geometry  statistics 


Property 

Value 

Space  dimension 

3 

Number  of  domains 

111 

Number  of  boundaries 

666 

Number  of  edges 

1472 

Number  of  vertices 

958 
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3-D  Fabric  Geometry 


2.2.1.  Warp  Spine  (wpl) 
Plane  Geometry  (wpl) 


Plane  Geometry’ 

Units 


Length  unit 

m 

Angular  unit 

deg 

Geometry  statistics 


Property 

Value 

Space  dimension 

2 

Number  of  domains 

0 

Number  of  boundaries 

0 

Center  Line  (pel) 

Position 
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3-D  Fabric  Geometry 


Name 

Value 

Position 

{0,0} 

Rotation  angle 


Name 

Value 

Maximum 

nx_fib*R_fib*2*pi 

Expressions 

{s,  R  fib * sin(s/(2  *R_fib)) } 

2.2.2.  Warp  Face  (wp2) 
Selections  of  resulting  entities 


Name 

Value 

Plane  type 

General 

Points 

{{0,  0,0},  {0,  1,0},  {-0.5,0,  1}} 

Draw  on  work  plane  in  3D 

On 

Plane  Geometry  (wp2) 


Plane  Geometry 
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3-D  Fabric  Geometry 


Units 


Length  unit 

m 

Angular  unit 

deg 

Geometry  statistics 


Property 

Value 

Space  dimension 

2 

Number  of  domains 

0 

Number  of  boundaries 

0 

Circle  1  (cl) 

Position 

Name 

Value 

Position 

{0,0} 

Radius 

R  fib 

2.2.3.  Weft  Spine  (wp3) 

Selections  of  resulting  entities 

Name 

Value 

Draw  on  work  plane  in  3D 

On 
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3-D  Fabric  Geometry 


Plane  Geometry  (wp3) 


Plane  Geometry 

Units 


Length  unit 

m 

Angular  unit 

deg 

Geometry  statistics 


Property 

Value 

Space  dimension 

2 

Number  of  domains 

0 

Number  of  boundaries 

0 

Center  Line  (pel) 

Position 


Name 

Value 

Position 

{0,0} 
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3-D  Fabric  Geometry 


Rotation  angle 


Name 

Value 

Maximum 

ny_fib*R_fib*2*pi 

Expressions 

{s,  -R_fib * sin(s/(2 *R_fib)) } 

2.2.4.  Weft  Face  (wp4) 
Selections  of  resulting  entities 


Name 

Value 

Plane  type 

General 

Points 

{{0,  0,  0},  {0,  1,0},  {0.5,  0,  1}} 

Draw  on  work  plane  in  3D 

On 

Plane  Geometry  (wp4) 


Plane  Geometry 


Units 


Length  unit 


m 
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3-D  Fabric  Geometry 


Angular  unit 


deg 


Geometry  statistics 


Property 

Value 

Space  dimension 

2 

Number  of  domains 

0 

Number  of  boundaries 

0 

Circle  1  (cl) 

Position 


Name 

Value 

Position 

{0,  0} 

Radius 

Rfib 

2.2.5.  Warp  Fiber  1  (swel) 
Selections  of  resulting  entities 


Name 

Value 

Keep  input  objects 

Off 

Include  all  inputs  in  finalize  operation 

Off 

2.2.6.  Weft  Fiber  1  (swe2) 

Selections  of  resulting  entities 

Name 

Value 

Keep  input  objects 

Off 

Include  all  inputs  in  finalize  operation 

Off 

Face-spine  alignment 

Adjust  spine 

2.2.7.  Rotate  Weft  Fiber  1  (rot2) 
Selections  of  resulting  entities 
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3-D  Fabric  Geometry 


Name 

Value 

Rotation 

90 

Point  on  axis  of  rotation 

{2*R_fib*pi/2,  0,  0} 

Axis  type 

z  -  axis 

2.2.8.  Warp  Fiber  2  (copyl) 
Selections  of  resulting  entities 


Name 

Value 

Keep  input  objects 

On 

X 

0 

y 

12.566370614359172 

z 

0 

2.2.9.  Flip  Warp  Fiber  2  (rotl) 
Selections  of  resulting  entities 


Name 

Value 

Rotation 

180 

Point  on  axis  of  rotation 

{0,  2*R_fib*pi,  0} 

Axis  type 

x  -  axis 

2.2.10.  Weft  Fiber  2  (copy3) 
Selections  of  resulting  entities 


Name 

Value 

Keep  input  objects 

On 

X 

12.566370614359172 

y 

0 

z 

0 
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3-D  Fabric  Geometry 


2.2.11.  Flip  Weft  Fiber  1  (rot3) 
Selections  of  resulting  entities 


Name 

Value 

Rotation 

180 

Point  on  axis  of  rotation 

{3*R_fib*pi,  0,  0} 

Axis  type 

y  -  axis 

2.2.12.  Warp(arrl) 
Selections  of  resulting  entities 


Name 

Value 

Array  type 

Linear 

Size 

ny_fib/2 

Displacement 

{0,  4*R_fib*pi,  0} 

2.2.13.  Weft  (arr2) 

Selections  of  resulting  entities 


Name 

Value 

Array  type 

Linear 

Size 

nx  fib/2 

Displacement 

{4*R_fib*pi,  0,  0} 

2.2.14.  Fabric  Layers  (arr3) 
Selections  of  resulting  entities 


Name 

Value 

Array  type 

Linear 

Size 

n  layer 

Displacement 

{0,  0,  4*R_fib} 
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3-D  Fabric  Geometry 


2.2.15.  Polymer  Matrix  (blkl) 
Position 


Name 

Value 

Position 

{0  -  (2*R  fib),  -(R_fib*pi)  -  (2*R_fib),  -(2*R_fib)  -  (2*R_fib)} 

Size  and  shape 


Name 

Size 


Value 

{2*pi*nx_fib*R_fib  +  (4*R_fib),  2*pi*ny_fib*R_fib  +  (4*R_fib),  4*R_fib*n_layer  + 
(4*R  fib)! 
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Appendix  D:  COMSOL  Model  Reports 


There  are  four  COMSOL-generated  model  reports  given.  One  for  the  Steam/ 
Propane  Reactor  Tube  model  and  one  each  for  the  three  parts  of  the  SiC  CVD  model: 
Flow  Model,  Chemical  Reaction  Model,  and  Moving  Mesh  Model. 
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Steam/Propane  Reactor  Tube 


1.  Global  Definitions 
1.1.  Parameters  1 
Parameters 


Name 

Expression 

Description 

R  fib 

.002 

fiber  radius  [m] 

N_fib 

10 

number  of  fibers  per  layer 

L_fib 

4 

number  of  layers 

u_in 

.05[m/s] 

inflow  velocity  [m/s] 

pin 

75  [Pa] 

Pressure  difference 

p_ref 

le5[Pa] 

Reference  pressure 

T_in 

700[K] 

Inlet  temperature 

ht 

100[W/mA2/K] 

Heat  transfer  coefficient 

DH2C3H8 

2.7e-5[mA2/s] 

Binary  diffusion  coefficient 

DC02C3H8 

5.1e-6[mA2/s] 

Binary  diffusion  coefficient 

DC3H8H20 

8.4e-6[mA2/s] 

Binary  diffusion  coefficient 

DC02H2 

3.6e-5[mA2/s] 

Binary  diffusion  coefficient 

DH2H20 

4.9e-5[mA2/s] 

Binary  diffusion  coefficient 

DC02H20 

l.le-5[mA2/s] 

Binary  diffusion  coefficient 

MC3H8 

44.1e-3[kg/mol] 

Molar  mass 

M_H2 

2.016e-3[kg/mol] 

Molar  mass 

MC02 

44.01e-3[kg/mol] 

Molar  mass 

MH20 

18.016e-3[kg/mol] 

Molar  mass 

w_C3H8_in 

0.28[1] 

Initial  weight  fraction 

w  H2  in 

0.01[1] 

Initial  weight  fraction 

w_C02_in 

0.01[1] 

Initial  weight  fraction 

H_sr 

410e3[J/mol] 

Enthalpy  of  reaction 

Cp 

2800[J/kg/K] 

Heat  capacity,  fluid 
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Steam/Propane  Reactor  Tube 


2.  Model  1  (modi) 

2.1.  Definitions 

2.1.1.  Variables 
Variables 
Selection 


Geometric  entity  level 


Entire  model 


Name  Expression 

Description 

Tout 

aveopl(T2) 

Average  outlet 
temperature 

rate 

chcs.x_w_C3H8*p/T2/R_const*7e5[l/s]*exp(-83.14e3[J/ 
mol]/R  const/T2) 

Reaction  rate 

dens_sr 

p/T2/R  const*(chcs.x  w  C3H8*M  C3H8  + 
chcs.x  w  H2*M  H2  +  chcs.x  w  C02*M  C02  + 
chcs.x  w  H20*M  H20) 

Density,  fluid 

2.1.2.  Model  Couplings 

Average  1 


Coupling  type 

Average 

Operator  name 

aveop 1 

Source  selection 


Geometric  entity  level 

Boundary 

Selection 

Boundary  4 

2.1.3.  Coordinate  Systems 
Boundary  System  1 


Coordinate  system  type 

Boundary  system 

Identifier 

sysl 
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Steam/Propane  Reactor  Tube 


Settings 


Name 

Value 

Coordinate  names 

{tl,  n,  to} 

Create  first  tangent  direction  from 

Global  Cartesian 

Length  unit 

m 

Angular  unit 

deg 

Geometry  statistics 


Property 

Value 

Space  dimension 

2 

Number  of  domains 

36 

Number  of  boundaries 

144 

2.2.  Geometry  1 
003$ 

0  03 
0.02S 
0  02 
0  01S 
0  01 
0  005 
0 

-COOS' 

•0  01 
•OOlS 

0  0  01  0  02  0  03  0  04  0  0) 

Geometry  1 

Units 


°o0o0o0o°o 

ououououou 
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Steam/Propane  Reactor  Tube 


2.2.1.  Circle  1  (cl) 
Position 


Name 

Value 

Position 

{0,0} 

Radius 

R  fib 

2.2.2.  Array  1  (arrl) 
Selections  of  resulting  entities 


Name 

Value 

Create  selections 

On 

Size 

{N_fib/2,  L  fib} 

Displacement 

{6*R_fib,  3*R  fib} 

2.2.3.  Copy  1  (copyl) 
Selections  of  resulting  entities 


Name 

Value 

Keep  input  objects 

On 

X 

0.0060 

y 

0.0030 

2.2.4.  Array  2  (arr2) 
Selections  of  resulting  entities 


Name 

Value 

Size 

{N_fib/2,  L_fib  -  1} 

Displacement 

{6*R_fib,  3*R  fib} 

2.2.5.  Rectangle  1  (rl) 
Position 
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Steam/Propane  Reactor  Tube 


Name 

Value 

Position 

{-2*R_fib,  -2*R_fib} 

Width 

N_fib*3*R_fib  +  Rfib 

Height 

3*R_fib*L  fib  +  Rfib 

Size 

{N_fib*3*R_fib  +  R  fib,  3*R_fib*L_fib  +  R_fib} 

2.3.  Materials 

2.3.1.  Pyrolytic  graphite  |soiid,as  deposited] 

00J 
0.025 
0.02 
0015 
0.01 
0  004 
0 

-0  004 
•0  01 

0  0  01  0  02  0  0)  004  0  04 


Pyrolytic  graphite  [solid, as  deposited] 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domains  2-36 
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Steam/Propane  Reactor  Tube 


2.4.  Laminar  Flow  (spf) 


Laminar  Flow 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Discretization  of  fluids 

PI  +P1 

Value  type  when  using  splitting  of 
complex  variables 

{Real,  Real,  Real,  Real,  Real,  Real,  Real,  Real,  Real} 

Equation  form 

Study  controlled 

Compressibility 

Compressible  flow  (Ma<0.3) 

Neglect  inertial  term  (Stokes  flow) 

Off 

Use  shallow  channel  approximation 

Off 

Allow  out  of  plane  properties 

1 
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Steam/Propane  Reactor  Tube 


Description 

Value 

Allow  radiation  properties 

1 

Allow  stokes  properties 

1 

Allow  turbulence  properties 

1 

Allow  MaxwellStefan  diffusion 

1 

Use  pseudo  time  stepping  for 
stationary  equation  form 

0 

Local  CFL  number 

1.3Amin(niterCMP,  9)  +  if(niterCMP>=25, 
9*1.3Amin(niterCMP  -  25,  9),  0)  +  if(niterCMP>=45, 
90*1.3Amin(niterCMP  -  45,  9),  0) 

Streamline  diffusion 

1 

Crosswind  diffusion 

1 

Isotropic  diffusion 

0 

Turbulence  model  type 

None 

Smoothing  parameter 

0.1 

_ 

Show  equation  assuming 

stdl/stat 
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Steam/Propane  Reactor  Tube 


2.4.1.  Fluid  Properties  1 


Fluid  Properties  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Density 

User  defined 

Density 

1 

Dynamic  viscosity 

User  defined 

Dynamic  viscosity 

1 

Reference  length 

1 

Reference  length  scale 

Automatic 

Mixing  length  limit 

Automatic 
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Steam/Propane  Reactor  Tube 


2.4.2.  Wall  1 


0 

02S 

0  0? 

01S 

o  or 

OOT 

0 

ooo 

oooo 

ooo 

oooo 

ooo 

oooo 

ooo 

oooo 

ooo 

oooo 

ooS 

0  01 

0  0  01  002  0  0)  0  04  0  05 

Wall  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundaries  2-3,  5-144 

Settings 


Description 

Value 

Temperature 

User  defined 

Temperature 

293. 1 5  [K] 

Electric  field 

User  defined 

Electric  field 

{0,  0,  0} 

Boundary  condition 

No  slip 

Apply  reaction  terms  on 

Individual  dependent  variables 

Use  weak  constraints 

0 
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Steam/Propane  Reactor  Tube 


2.4.3.  Initial  Values  1 


0  0  01  002  0.0S  004  0  05 

Initial  Values  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Velocity  field 

{0,  0,  0} 

Pressure 

pref 

0  03 
e  o iS 


u 

0  01 
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2.4.4.  Inlet  1 


Inlet  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundary  1 

Settings 


Description 

Value 

Apply  reaction  terms  on 

All  physics  (symmetric) 

Use  weak  constraints 

0 

Boundary  condition 

Pressure,  no  viscous  stress 

Pressure 

pin  +  pref 

Standard  pressure 

1  [atm] 

Standard  molar  volume 

0.0224 136[mA3/mol] 

Normal  mass  flow  rate 

le-5[kg/s] 
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Description 

Value 

Mass  flow  type 

Mass  flow  rate 

Standard  flow  rate  defined  by 

Standard  density 

Channel  thickness 

1.0  [m] 

2.4.5.  Outlet  1 


Outlet  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundary  4 

Settings 


Description 

Value 

Apply  reaction  terms  on 

All  physics  (symmetric) 

Use  weak  constraints 

0 

Boundary  condition 

Pressure,  no  viscous  stress 

Pressure 

pref 
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2.5.  Heat  Transfer  in  Fluids  2  (ht2) 

0  03 
0.02 £ 

0  02 
0  015 
0  Of 
0005 
0 

-oooS 
■0  01 

0  0  01  002  0  03  0  04  0  05 

Heat  Transfer  in  Fluids  2 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Temperature 

Linear 

Compute  boundary  fluxes 

1 

Apply  smoothing  to  boundary  fluxes 

1 

Value  type  when  using  splitting  of  complex  variables 

Real 

Equation  form 

Study  controlled 

Out-of-plane  heat  transfer 

0 

Surface-to-surface  radiation 

0 

Radiation  in  participating  media 

0 
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Description 

Value 

Heat  transfer  in  biological  tissue 

0 

Heat  transfer  in  porous  media 

0 

Moving  frame 

Moving  frame  flag  (false) 

Streamline  diffusion 

1 

Crosswind  diffusion 

1 

Isotropic  diffusion 

0 

Show  equation  assuming 

stdl/stat 

2.5.1.  Heat  Transfer  in  Fluids  1 


0  0  01  002  0  03  0  04  0  03 

Heat  Transfer  in  Fluids  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 
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Settings 


Description 

Value 

Fluid  type 

Gas/Liquid 

Thermal  conductivity 

User  defined 

Thermal  conductivity 

{{k_sr,  0,  0},  {0,  k_sr,  0},  {0,  0,  k_sr} } 

Density 

User  defined 

Density 

1 

Fleat  capacity  at  constant  pressure 

User  defined 

Fleat  capacity  at  constant  pressure 

Cp 

Ratio  of  specific  heats 

User  defined 

Ratio  of  specific  heats 

1.5 

2.5.2.  Thermal  Insulation  1 


Thermal  Insulation  1 


Selection 


Geometric  entity  level 


Boundary 
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Selection 


Boundary  2 


2.5.3.  Initial  Values  1 


Initial  Values  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Temperature 

293. 1 5  [K] 
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2.5.4.  Surface-to-Ambient  Radiation  1 


Surface-to- Ambient  Radiation  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundary  3 

Settings 


Description 

Value 

Ambient  temperature 

293. 1 5  [K] 

Surface  emissivity 

User  defined 

Surface  emissivity 

.998 
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2.5.5.  Outflow  1 

0  02f 
0  0? 

0  018 
0  018 
0  014 
0  012 
0  01 
0  ootfj 
0  0M 
0  004" 

0  00? 

a 

0  00? 

0  004^ 

004  0  045  00)  00))  000  0  00)  007 


Outflow  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundary  4 

00 
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2.5.6.  Boundary  Heat  Source  1 


Boundary  Heat  Source  1 


Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundary  2 

Settings 


Description 

Value 

Heat  source 

General  source 

Boundary  heat  source 

User  defined 

Boundary  heat  source 

1000 
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2.5.7.  Heat  Flux  1 


Heat  Flux  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundaries  5-144 

Settings 


Description 

Value 

Heat  flux 

General  inward  heat  flux 

Inward  heat  flux 

ht*(Tl  -T2) 

Heat  transfer  coefficient 

-5 

External  temperature 

T1 
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2.5.8.  Heat  Source  1 


Heat  Source  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Heat  source 

General  source 

Heat  source 

User  defined 

Heat  source 

-rate*H  sr 
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2.5.9.  Temperature  1 


Temperature  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundary  1 

Settings 


Description 

Value 

Temperature 

600 

Classic  constraints 

Apply  reaction  terms  on 

All  physics  (symmetric) 

Use  weak  constraints 

0 
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2.6.  Heat  Transfer  in  Solids  (ht) 

0  03 

0.026 
0  02 
0016 
001 
0  006 
0 

-0  006 
001 

u  0  0  006  0  01  0  016  0  02  0  025  0  03  0  035  0  04  0  045  0  05  C 

Heat  Transfer  in  Solids 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domains  2-36 

Heat  Transfer  in  Solids  (ht) 


Settings 


Description 

Value 

Temperature 

Quadratic 

Compute  boundary  fluxes 

1 

Apply  smoothing  to  boundary  fluxes 

1 

Value  type  when  using  splitting  of  complex  variables 

Real 

Equation  form 

Study  controlled 

Out-of-plane  heat  transfer 

0 

Surface-to-surface  radiation 

0 
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Description 

Value 

Radiation  in  participating  media 

0 

Heat  transfer  in  biological  tissue 

0 

Heat  transfer  in  porous  media 

0 

Moving  frame 

Moving  frame  flag  (false) 

Streamline  diffusion 

1 

Crosswind  diffusion 

1 

Isotropic  diffusion 

0 

Show  equation  assuming 

stdl/stat 

2.6.1.  Heat  Transfer  in  Solids  1 


Heat  Transfer  in  Solids  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domains  2-36 
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Settings 

Description 

Value 

Thermal  conductivity 

From  material 

Thermal  conductivity 

{{0,  0,  0},  {0,  0,0},  {0,  0,  0}} 

Density 

From  material 

Heat  capacity  at  constant  pressure 

From  material 

Properties  from  material 


Property 

Material 

Property  group 

Thermal  conductivity 

Pyrolytic  graphite  [sol id, as  deposited] 

Basic 

Density 

Pyrolytic  graphite  [sol id, as  deposited] 

Basic 

Heat  capacity  at  constant 

Pyrolytic  graphite  [sol id, as  deposited] 

Basic 

pressure 

2.6.2.  Thermal  Insulation  1 
Selection 


Geometric  entity  level 

Boundary 

Selection 

No  boundaries 
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2.6.3. 

0  03 

0.02S 
0  02 
001S 
0  01 
0  006 
0 

-OOOS 
001 

0  0  005  0  01  0  015  0  02  0  025  0  01  0  015  0  04  0  045  0  05  C 

Initial  Values  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domains  2-36 

Initial  Values  1 


Settings 


Description 

Value 

Temperature 

293. 1 5  [K] 

87 


Steam/Propane  Reactor  Tube 


2.6.4.  Heat  Flux  1 


Heat  Flux  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundaries  5-144 

Settings 


Description 

Value 

Heat  flux 

General  inward  heat  flux 

Inward  heat  flux 

-ht*(Tl  -T2) 

Heat  transfer  coefficient 

5 

External  temperature 

T2 
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2.7.  Transport  of  Concentrated  Species  (chcs) 


Transport  of  Concentrated  Species 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Mass  fraction 

Linear 

Compute  boundary  fluxes 

1 

Apply  smoothing  to  boundary  fluxes 

1 

Value  type  when  using  splitting  of  complex  variables 

Real 

Equation  form 

Study  controlled 

Migration  in  electric  field 

0 

Convection 

1 
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Description 

Value 

Diffusion  model 

Mixture-averaged 

Residual 

Approximate  residual 

From  mass  constraint 

wH20 

Streamline  diffusion 

1 

Crosswind  diffusion 

1 

Crosswind  diffusion  type 

Do  Carmo  and  Galeao 

Isotropic  diffusion 

0 

Regularization 

On 

Enable  space-dependent  physics  interfaces 

0 

Synchronize  with  COMSOL  Multiphysics 

Show  equation  assuming 

stdl/stat 

2.7.1.  Convection  and  Diffusion 


Convection  and  Diffusion 
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Selection 


1 

Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Velocity  field 

Velocity  field  (spf/fpl) 

Electric  potential 

User  defined 

Electric  potential 

0 

Thermal  diffusion 
coefficient 

{0, 0,  0,  0} 

Molar  mass 

{M  H20,  M  C3H8,  M_H2,  M  C02} 

Density 

denssr 

Mixture  density 

User  defined 

Maxwell-Stefan  diffusivity 
matrix 

{{1,  D  C3H8  H20,  D  H2  H20,  D  C02  H20}, 

{D  C3H8  H20,  1,  D  H2  C3H8,  D  C02  C3H8},  {D  H2  H20, 
D  H2  C3H8,  1,  D  C02  H2},  {D  C02  H20,  D  C02  C3H8, 
DC02H2,  1}} 

Diffusion  coefficient 

{{{le-5[mA2/s],  0,  0},  {0,  le-5[mA2/s],  0},  {0,  0,  le-5[mA2/s]}}, 
{{le-5[mA2/s],  0,  0},  {0,  le-5[mA2/s],  0},  {0,  0,  le-5[mA2/s]}}, 
{{le-5[mA2/s],  0,  0},  {0,  le-5[mA2/s],  0},  {0,  0,  le-5[mA2/s]}}, 
{{le-5[mA2/s],  0,  0},  {0,  le-5[mA2/s],  0},  {0,  0,  le-5[mA2/s]}}} 
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2.7.2.  No  Flux  1 


0 

02* 

0  0? 

019 

0  of 

009 

0 

ooo 

oooo 

ooo 

oooo 

ooo 

oooo 

ooo 

oooo 

ooo 

oooo 

00* 

0  01 

0  0  01  002  0  0)  0  04  0  09 


Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundaries  2-3,  5-144 

Settings 


Description 

Value 

Apply  for  all  species 

Apply  for  all  species 
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2.7.3.  Initial  Values  1 


Initial  Values  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Mass  fraction 

w_C3H8_in 

Mass  fraction 

1/2 

Mass  fraction 

w  H2  in 

Mass  fraction 

w_C02_in 
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2.7.4.  Reactions  1 


Reactions  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Total  rate  expression 

{0,  -M  C3H8*rate,  10*M  H2*rate,  3*M  C02*rate} 

Mass  transfer  to  other  phases 

0 
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2.7.5.  Inflow  1 


Inflow  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundary  1 

Settings 


Description 

Value 

Mixture  specification 

Mass  fractions 

Mass  fraction 

{0,  w C3H8 in,  w H2 in,  w C02 in} 

Apply  reaction  terms  on 

All  physics  (symmetric) 

Use  weak  constraints 

0 
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2.7.6.  Outflow  1 

0  0221 
0  0? 

0  018 
0  018 
0  014 
0  012 
0  01 
OOOtfj 

0  0M 
0  004“ 

0  00? 

d 

000? 

•0  004 

004  0  045  00)  00))  000  0  00)  007 


()  () 
'  o  ^ 

0X0 

o 

Cj  Cj 
(  ) 


Outflow  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundary  4 

2.8.  Mesh  1 
Mesh  statistics 


Property 

Value 

Minimum  element  quality 

0.6495 

Average  element  quality 

0.9241 

Triangular  elements 

6484 

Edge  elements 

666 

Vertex  elements 

144 

00 
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Mesh  1 


2.8.1.  Size  (size) 
Settings 


Name 

Value 

Maximum  element  size 

0.0062 

Minimum  element  size 

1.24E-4 

Resolution  of  curvature 

0.4 

Maximum  element  growth  rate 

1.4 

Predefined  size 

Coarse 

2.8.2.  Free  Triangular  3  (ftri3) 
Selection 


Geometric  entity  level 


Remaining 
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3.  Study  1 
3.1.  Stationary 
Study  settings 


Property 

Value 

Include  geometric  nonlinearity 

Off 

Mesh  selection 


Geometry 

Mesh 

Geometry  1  (geoml) 

meshl 

Physics  selection 


Physics  Discretization 

Laminar  Flow  (spf) 

physics 

Heat  Transfer  in  Fluids  2  (ht2) 

physics 

Heat  Transfer  in  Solids  (ht) 

physics 

Transport  of  Concentrated  Species  (chcs) 

physics 

3.2.  Solver  Configurations 
3.2.1.  Solver  1 

Compile  Equations:  Stationary  (stl) 
Study  and  step 


Name 

Value 

Use  study 

Study  1 

Use  study  step 

Stationary 

Dependent  Variables  1  (vl) 
General 


Name 

Value 

Defined  by  study  step 

Stationary 
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Initial  values  of  variables  solved  for 


Name 

Value 

Solution 

Zero 

Values  of  variables  not  solved  for 


Name 

Value 

Solution 

Zero 

Temperature  (modl.T2)  (modl_T2) 

General 


Name 

Value 

Field  components 

modl.T2 

Velocity  field  (modl.u)  (modl  u) 

General 


Name 

Value 

Field  components 

{modl.u,  modl.v} 

Temperature  (modl.Tl)  (modl_Tl) 

General 


Name 

Value 

Field  components 

modl.Tl 

Mass  fraction  (modl.w_C3H8)  (modl_w_C3H8) 

General 


Name 

Value 

Field  components 

modl.w_C3H8 

Pressure  (modl.p)  (modl_p) 

General 
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Name 

Value 

Field  components 

modl.p 

Mass  fraction  (modl.w_C02)  (modl_w_C02) 

General 


Name 

Value 

Field  components 

modl.w_C02 

Mass  fraction  (modl.w_H2)  (modl_w_H2) 

General 


Name 

Value 

Field  components 

modl.w  F12 

Stationary  Solver  1  (si) 
General 


Name 

Value 

Defined  by  study  step 

Stationary 

Fully  Coupled  1  (fcl) 


General 


Name 

Value 

Linear  solver 

Direct  1 

Direct  1  (dl) 

General 


Name 

Value 

Solver 

PARD1SO 
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4.  Results 

4.1.  Data  Sets 

4.1.1.  Solution  1 
Selection 


Geometric  entity  level 

Domain 

Selection 

Geometry  geoml 

Solution 


Name 

Value 

Solution 

Solver  1 

Model 

Save  Point  Geometry  1 

4.2.  Plot  Groups 
4.2.1.  Velocity  (spf) 

Sorfaco  Velocity  magnitude  (rrVsl  Streamline  Velocity  field 


0  03 
0  035 
0  02 
0  015 
0  01 
0  00b 
0 

•0  005 
0.01 
•0  015 


0  01 


0  02 


0  02 


004 


005 


▲  1  7744K10’3 

1  0  0016 
0  0014 
0  0012 
0  001 
00008 
00006 

100004 
0  0002 
0 


Surface:  Velocity  magnitude  (m/s)  Streamline:  Velocity  field 
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4.2.2.  Pressure  (spf) 

Contour:  Pressure  iPa)  Contour  Pressure  (Pa) 


0  045 
004 
0  015 
0  03 
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0.015 
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Contour:  Pressure  (Pa)  Contour:  Pressure  (Pa) 


4.2.3.  Temperature 

Surface  Tomperature  00  Surface  Temperature  HO 


005 
0  045 
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4.2.4.  Mass  Fractions 


Lin#  Graph:  Mats  fraction  Mass  fraction  (1)  lata  Graph  Mass  fraction  (1) 
Lin*  Graph  Mass  fraction  (1)  line  Graph  Mass  fraction  (1> 


Line  Graph:  Mass  fraction  Mass  fraction  (1)  Line  Graph:  Mass  fraction  (1)  Line  Graph:  Mass  fraction  (1) 
Line  Graph:  Mass  fraction  (1) 
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1.  Global  Definitions 
1.1.  Parameters  1 
Parameters 


Name 

Expression  Description 

R  fib 

.002 

fiber  radius  [m] 

N_fib 

6 

number  of  fibers  per  layer 

L_fib 

4 

number  of  layers 

u_in 

.05[m/s] 

inflow  velocity  [m/s] 

p_in 

3  [Pa] 

Pressure  difference 

p_ref 

le5[Pa] 

Reference  pressure 

T_in 

700[K] 

Inlet  temperature,  reformer  bed 

ht 

100[W/mA2/K] 

Heat  transfer  coefficient,  heating  tubes 

DH2C3H8 

2.7e-5[mA2/s] 

Binary  diffusion  coefficient 

DC02C3H8 

5.1e-6[mA2/s] 

Binary  diffusion  coefficient 

DC3H8H20 

8.4e-6[mA2/s] 

Binary  diffusion  coefficient 

DC02H2 

3.6e-5[mA2/s] 

Binary  diffusion  coefficient 

DH2H20 

4.9e-5[mA2/s] 

Binary  diffusion  coefficient 

DC02H20 

l.le-5[mA2/s] 

Binary  diffusion  coefficient 

MC3H8 

44.1e-3[kg/mol] 

Molar  mass 

M_H2 

2.016e-3[kg/mol] 

Molar  mass 

MC02 

44.01e-3[kg/mol] 

Molar  mass 

MH20 

1 8.0 16e-3  [kg/mol] 

Molar  mass 

w_C3H8_in 

0.28[1] 

Initial  weight  fraction 

w  H2  in 

0.01[1] 

Initial  weight  fraction 

w_C02_in 

0.01[1] 

Initial  weight  fraction 

H_sr 

410e3[J/mol] 

Enthalpy  of  reaction 

Cp 

2800[J/kg/K] 

Heat  capacity,  fluid 

k  sr 

0.1[W/m/K] 

Thermal  conductivity,  reformer  bed 
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Name 

Expression 

Description 

visc_sr 

2.7e-5[Pa*s] 

Viscosity,  reformer  bed 

X_disp 

3*R_fib 

displacement  in  x  direction  between  fiber  centers 

Y_disp 

3*R_fib 

displacement  in  y  direction  between  fiber  centers 
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2.  Model  1  (modi) 

2.1.  Definitions 

2.1.1.  Variables 
Variables 
Selection 

Geometric  entity  level  Entire  model 


Name 

Expression 

Description 

Tout 

aveop(T2) 

Average  outlet 
temperature 

rate 

chcs.x_w_C3H8*p/T2/R_const*7e5[l/s]*exp(-83.14e3[J/mol]/ 

R  const/T2) 

Reaction  rate 

denssr 

p/T2/R  const*(chcs.x  w  C3H8*M  C3H8  +  chcs.x  w  H2*M  H2  + 
chcs.x_w_C02*M_C02  +  chcs.x_w_H20*M_H20) 

Density,  reformer  bed 

U_out 

aveop(u) 

Average  outlet  velocity 

Rema 

X 

maxop  1  (spf.cellRe) 

Maximum  Cell 

Reynolds  number 

2.1.2.  Model  Couplings 

Average 


Coupling  type 

Average 

Operator  name 

aveop 

Source  selection 


Geometric  entity  level 

Boundary 

Selection 

Boundary  7 

Maximum  1 


Coupling  type 

Maximum 

Operator  name 

maxop  1 
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Source  selection 


Geometric  entity  level 

Domain 

Selection 

Domain  2 

2.1.3.  Selections 
Carbon  Fibers 

Selection  type 

Explicit 

Selection 

Domains  3-23 


0  005  001 


0015  0  02  0  025  0  0 


Carbon  Fibers 
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Fluid  Matrix 


Selection  type 

Explicit 


Selection 


No  domains 


Inlet 


Selection  type 

Explicit 


Selection 


Boundary  7 


Inlet 
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Outlet 

Selection  type 

Explicit 

1 

Selection 

No  boundaries 


Air  Boundary 

Selection  type 

Explicit 

Selection 

No  boundaries 


Hot  Plate 

Selection  type 

Explicit 

1 - 

Selection 

No  boundaries 


2.1.4.  Coordinate  Systems 
Boundary  System  1 


Coordinate  system  type 

Boundary  system 

Identifier 

sysl 

Settings 
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Name 

Value 

Coordinate  names 

{tl,  n,  to} 

Create  first  tangent  direction  from 

Global  Cartesian 

2.2.  Geometry  1 


Geometry’  1 

Units 


Length  unit 

m 

Angular  unit 

deg 

Geometry  statistics 


Property 

Value 

Space  dimension 

2 

Number  of  domains 

23 

Number  of  boundaries 

91 
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2.2.1.  Circle  1  (cl) 
Position 


Name 

Value 

Position 

{0,0} 

Radius 

Rfib 

2.2.2.  Array  1  (arrl) 
Selections  of  resulting  entities 


Name 

Value 

Create  selections 

On 

Size 

{N_fib/2,  Lfib} 

Displacement 

{2*X_disp,  Y  dispj 

2.2.3.  Copy  1  (copyl) 

Selections  of  resulting  entities 

Name 

Value 

Keep  input  objects 

On 

X 

0.0060 

y 

0.0030 

2.2.4.  Array  2  (arr2) 
Selections  of  resulting  entities 


Name 

Value 

Size 

{N_fib/2,  Lfib  -  1 } 

Displacement 

{2*X_disp,  Y_disp} 

2.2.5.  Rectangle  1  (rl) 
Position 
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Name 

Value 

Position 

{-2*R_fib,  -2*R_fib} 

Width 

N_fib*X_disp  +  Rfib 

Height 

L_fib*Y_disp  +  Rfib 

Size 

{N_fib*X_disp  +  R  fib,  L_fib*Y_disp  +  R  fib} 

2.2.6.  Rectangle  2  (r2) 
Position 


Name 

Value 

Position 

{-.03,  -2*R_fib} 

Width 

(N_fib*X_disp  -  R  fib)  +  0.03 

Height 

L_fib*Y_disp  +  Rfib 

Size 

{(N_fib*X_disp  -  R  fib)  +  0.03,  L_fib*Y_disp  +  R  fib} 

2.3. 

2.3.1. 

0.02 
001B 
0016 
0.014 
001? 

0.01 
0  006 
0  006 
0  004 
000? 

0 

- 

0  0  005  001  0015  002  0  025  0  0 

Pyrolytic  graphite  [solid, as  deposited] 


Materials 

Pyrolytic  graphite  |  solid, as  deposited] 
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Selection 


Geometric  entity  level 

Domain 

Selection 

Domains  3-23 

2.4.  Laminar  Flow  (spf) 


Laminar  Flow 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domains  1-2 

Settings 


Description 

Value 

Discretization  of  fluids 

PI  +P1 

Value  type  when  using  splitting  of  complex 
variables 

{Real,  Real,  Real,  Real,  Real,  Real,  Real, 

Real,  Real} 

Equation  form 

Study  controlled 
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Description 

Value 

Compressibility 

Compressible  flow  (Ma<0.3) 

Neglect  inertial  term  (Stokes  flow) 

Off 

Use  shallow  channel  approximation 

Off 

Allow  out  of  plane  properties 

1 

Allow  radiation  properties 

1 

Allow  stokes  properties 

1 

Allow  turbulence  properties 

1 

Allow  MaxwellStefan  diffusion 

1 

Use  pseudo  time  stepping  for  stationary 
equation  form 

0 

Local  CFL  number 

1.3Amin(niterCMP,  9)  +  if(niterCMP>=25, 
9*1.3Amin(niterCMP  -  25,  9),  0)  + 
if(niterCMP>=45,  90*1.3Amin(niterCMP  - 
45,  9),  0) 

Streamline  diffusion 

1 

Crosswind  diffusion 

1 

Isotropic  diffusion 

0 

Turbulence  model  type 

None 

Smoothing  parameter 

0.1 

Show  equation  assuming 

stdl/stat 
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2.4.1.  Fluid  Properties  1 


•002 


OOl 


0  01 


002 


Fluid  Properties  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domains  1-2 

Settings 


Description 

Value 

Density 

User  defined 

Density 

1000 

Dynamic  viscosity 

User  defined 

Dynamic  viscosity 

8.9e-4 

Reference  length 

1 

Reference  length  scale 

Automatic 

Mixing  length  limit 

Automatic 
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2.4.2.  Wall  1 


0  0? 


0  023 


0  02 


0015 


0  01 


o  ooS 


•0  003 


0  01 


0015 


oooo 

ooo 

oooo 

ooo 

oooo 

ooo 

•0  02  -001 


0  01  002  0  05 


Wall  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundaries  2-3,  5-6,  8-91 

Settings 


Description 

Value 

Temperature 

User  defined 

Temperature 

293. 1 5  [K] 

Electric  field 

User  defined 

Electric  field 

{0,  0,  0} 

Boundary  condition 

No  slip 

Apply  reaction  terms  on 

Individual  dependent  variables 

Lise  weak  constraints 

0 
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2.4.3.  Initial  Values  1 


Initial  Values  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domains  1-2 

Settings 


Description 

Value 

Velocity  field 

{0,  0,  0} 

Pressure 

pref 
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2.4.4.  Inlet  1 


Inlet  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundary  1 

Settings 


Description 

Value 

Apply  reaction  terms  on 

All  physics  (symmetric) 

Use  weak  constraints 

0 

Boundary  condition 

Pressure,  no  viscous  stress 

Pressure 

p_in  +  p_ref 

Standard  pressure 

1  [atm] 

Standard  molar  volume 

0.0224 136[mA3/mol] 

Normal  mass  flow  rate 

le-5[kg/s] 
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Description 

Value 

Mass  flow  type 

Mass  flow  rate 

Standard  flow  rate  defined  by 

Standard  density 

Channel  thickness 

1.0[m] 

2.4.5.  Outlet  1 


Outlet  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundary  7 

Settings 


Description 

Value 

Apply  reaction  terms  on 

All  physics  (symmetric) 

Use  weak  constraints 

0 

Boundary  condition 

Pressure,  no  viscous  stress 
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Description 

Value 

Pressure 

pref 

2.5.  Mesh  1 
Mesh  statistics 


Property 

Value 

Minimum  element  quality 

0.1701 

Average  element  quality 

0.8061 

Triangular  elements 

8963 

Quadrilateral  elements 

1260 

Edge  elements 

716 

Vertex  elements 

90 

OOT 
0  02} 


■0  00} 
0.01 


0015* 


■0  02 


■0  01 


0  01 


002 


0  01 


Mesh  1 


2.5.1.  Size  (size) 
Settings 
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Name 

Value 

Maximum  element  size 

0.0064 

Minimum  element  size 

1.28E-4 

Resolution  of  curvature 

0.4 

Maximum  element 
growth  rate 

1.4 

Predefined  size 

Coarse 

2.5.2.  Size  1  (sizel) 
Selection 


Geometric  entity  level 

Domain 

Selection 

Domains  1-2 

0  0J 
0  02} 


-0  00 } 
0.01 


0  015  0  02  -0  01  0  0  01  002  OOS 

Size  1 

Settings 
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Name 

Value 

Calibrate  for 

Fluid  dynamics 

Maximum  element  size 

0.00226 

Minimum  element  size 

1.04E-4 

Resolution  of  curvature 

0.6 

Maximum  element  growth  rate 

1.25 

Predefined  size 

Coarser 

2.5.3.  Size  2  (size2) 
Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundaries  2-3,  5-6,  8-91 

4004 

0.01 

0015  0  02  0  01  0  0  01  0  02  0  01 


Size  2 

Settings 
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Name 

Value 

Calibrate  for 

Fluid  dynamics 

Maximum  element  size 

0.00117 

Minimum  element  size 

5.2E-5 

Resolution  of  curvature 

0.3 

Maximum  element  growth  rate 

1.15 

2.5.4.  Corner  Refinement  1  (crl) 
Selection 


0  02 
0  015 
0.01 
OOOff 
0 

■0  005 
0  01 
0015' 


1  • 

4  ►  4^  .4  k  -*></»■<» ; 

1158SistIP'; 


•0  02 


•OOl 


0  01 


002 


0  05 


Corner  Refinement  1 


2.5.5.  Free  Triangular  1  (ftril) 
Selection 


Geometric  entity  level 


Domain 
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Selection 


Domains  1-2 


Free  Triangular  1 

2.5.6.  Boundary  Layers  1  (bll) 
Selection 


Geometric  entity  level 

Domain 

Selection 

Domains  1-2 
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4  02 


4  01 


0  01 


002 


0  0) 


Boundary  Layers  1 

Settings 


Name 

Value 

Handling  of  shaip  comers 

Trimming 

Boundary  Layer  Properties  1  (blpl) 
Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundaries  2-3,  5-6,  8-91 
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4004 

0.01 

0015  0  02  0  01  0  0  01  0  02  0  0S 


Boundary  Layer  Properties  1 


Settings 


Name 

Value 

Number  of  boundary  layers 

2 

2.5.7.  Free  Triangular  2  (ftri2) 
Selection 


Geometric  entity  level 


Remaining 
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3.  Study  1 

3.1.  Parametric  Sweep 

Parameter  name:  p_in 

Parameters: 

3.2.  Stationary 
Study  settings 


Property 

Value 

Include  geometric  nonlinearity 

Off 

Mesh  selection 


Geometry 

Mesh 

Geometry  1  (geoml) 

meshl 

Physics  selection 

Physics 

Discretization 

Laminar  Flow  (spf) 

physics 

3.3.  Solver  Configurations 
3.3.1.  Solver  1 


Compile  Equations:  Stationary  (stl) 
Study  and  step 


Name 

Value 

Use  study 

Study  1 

Use  study  step 

Stationary 

Dependent  Variables  1  (vl) 
General 


Name 

Value 

Defined  by  study  step 

Stationary 
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Initial  values  of  variables  solved  for 


Name 

Value 

Solution 

Zero 

Values  of  variables  not  solved  for 


Name 

Value 

Solution 

Zero 

Velocity  field  (modl.u)  (modl  u) 

General 


Name 

Value 

Field  components 

{modl.u,  modl.v} 

Pressure  (modl.p)  (modl_p) 

General 

Name 

Value 

Field  components 

modl.p 

Stationary  Solver  1  (si) 
General 


Name 

Value 

Defined  by  study  step 

Stationary 

General 


Name 

Value 

Defined  by  study  step 

Parametric  Sweep 

Parameter  value  list 

range(0.01,  0.01,  0.5) 

_ 
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Fully  Coupled  1  (fcl) 


General 


Name 

Value 

Linear  solver 

Direct  1 

129 


Flow  Model 


4.  Results 

4.1.  Data  Sets 

4.1.1.  Solution  1 
Selection 


Geometric  entity  level 

Domain 

Selection 

Geometry  geoml 

Solution 


Name 

Value 

Solution 

Solver  1 

Model 

Save  Point  Geometry  1 

4.1.2.  Solution  2 
Selection 


Geometric  entity  level 

Domain 

Selection 

Geometry  geoml 

Solution 


Name 

Value 

Solution 

Parametric  2 

Model 

Save  Point  Geometry  1 

4.2.  Tables 

4.2.1.  Evaluation  2D 

Interactive  2D  values 

Evaluation  2D 


X 

y 

Value 

0.01189 

0.01756 

656.48334 
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4.3.  Plot  Groups 


4.3.1.  Velocity  (spf) 

pjntS0)-O.5  Surface  Velocity  magnitude  IrrVs)  Streamline  Velocity  field 

A  0  Olfif 

0  03 

“  0016 
0  014 
0  012 
001 
0  006 
0006 
0  004 
0  002 
0 


•0  02 


-0  01 


0  01 


0  02 


V  0 


p_in(50)=0.5  Surface:  Velocity  magnitude  (m/s)  Streamline:  Velocity  field 


4.3.2.  Pressure  (spf) 

p_in450)-0.5  Contour  Pressure  (Pa) 

0  035 
0  03 
0.025  - 
0  02 
0.015 
001 
0005 
0 

•0.005 
•0  01- 
•0  015 


a  i*io3 

XlO4 


-0  02 


•0  01 


001 


0  02  O03vixl05 


p_in(50)=0.5  Contour:  Pressure  (Pa) 
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4.3.3.  Outlet  Velocity 

Lin*  Graph  Averago  outlot  velocity  (rrVsl 


0  004 

0  0035 

*•» 

X  0  003 

§  0  0025* 

1 

•  0  002- 

i 

I  0  0015- 
<  0001 
0  0005- 


01  0  2  0  3  0  4  0  5 

Prew >jr*  <Mf««en<e  (Pa) 


Line  Graph:  Average  outlet  velocity  (m/s) 

4.3.4.  Max  Reynolds  Number 

Une  Graph:  Maximum  Call  Reynolds  number  (1) 


I 

“O 

? 


2  8 
2  b 
7  4 

7.7 

7 

18 

16- 
14  - 
12- 
1- 
06 
0  b 
0  4 
0.2  , 


01 


0  2  0  3 

Prew«xe  4*Oe*ence  (Pa) 


04 


05 


Line  Graph:  Maximum  Cell  Reynolds  number  (1) 
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4.3.5.  Velocity  (spf)  1 

Xj»np|58)-0  0194  Sur4»c#:  V#tocrty  magmtud*  I  m/s  I 

▲  0  047 


-001  -0  01 


0  03*  0 


X_disp(58)=0.0154  Surface:  Velocity  magnitude  (m/s) 


4.3.6.  Pressure  (spf)  1 

X_disp<98>-0  0154  Contour:  Pressure  |P») 


ooe 

0  09 

0  04  * 

0  02 

- 

002 

f',  i~\ 

K  °  n  °  K  0 

001 

K  O  o  0  0 

O 

R  o  g  °  g  o 

-0  01  - 

-1 

•0  02* 

* 

4)03 

•4  04 

•0  02  0  0  02  0  04  0  00  0  00 


▼  1x10* 


X_disp(58)=0.0154  Contour:  Pressure  (Pa) 
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1.  Global  Definitions 
1.1.  Parameters  1 
Parameters 


Name 

Expression 

Description 

R 

R  const 

L_r 

270[mm] 

length  of  reactor 

R_r 

135  [mm] 

radius  of  reactor 

V_r 

L_r*pi*R_rA2 

Volume  of  reactor 

A_r 

L_r*2*pi*R_r 

Surface  area  of  reactor 

P_r 

250[mbar] 

pressure  of  reactor 

T_o 

2073 [K] 

initial  temperature 

c_SiH4_o 

(3/12)*p_r/(R_const*T_o) 

inital  silane  concentration 

c  H2  o 

(8/12)*p_r/(R_const*T_o) 

initial  hydrogen  concentration 

c_C3H8_o 

(l/12)*p  r/(R  const*T  o) 

initial  propane  concnetration 

Temp 

2073 [K] 

Initial  temperature 
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2.  Model  1  (modi) 

2.1.  Definitions 

2.1.1.  Variables 


Variables  1 
Selection 


Geometric  entity  level 


Entire  model 


Name  Expression 

Description 

Si_gas 

modl.re.c  SiH4  +  modl.re.c  SiH2  + 
2*modl.re.c  Si2H6  +  modl.re.c  Si  + 
2*modl.re.c  Si2  +  2*modl.re.c  Si2C  + 
modl.re.c_SiCH2 

Total  Concentration  of  Si  in  Gas 
Form 

C_gas 

3*modl.re.c  C3H8  +  modl.re.c  CH3  + 
2*modl.re.c  C2H5  +  modl.re.c  CH4 + 
2*modl.re.c  C2H6  +  2*modl.re.c  C2H4  + 
2*modl.re.c  C2H2  +  modl.re.c  Si2C  + 
modl.re.c  SiCH2  +  modl.re.c  CH2 

Total  Concentration  of  C  in  Gas 
phase 

Si_solid 

c_SiH4_o  -  Si  gas 

Concentration  of  Si  atoms 
deposited  onto  surface 

Csolid 

3*c_C3H8_o  -  C  gas 

Concentration  of  C  atoms 
deposited  onto  surface 

2.2.  Reaction  Engineering  (re) 
Selection 


Geometric  entity  level 

Entire  model 

Settings 

Description 

Value 

Equation  form 

Study  controlled 
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Description 

Value 

Old  pressure 

R  const*T*(c  SiH4  +  c  SiH2  +  c  H2  +  c  Si2H6  +  c  Si  +  c  H  + 
c_C3H8  +  c_CH3  +  c_C2H5  +  c_CH4  +  c_C2H6  +  c_C2H4  + 
c  C2H2  +  c  Si2  +  c  Si2C  +  c  SiCH2  +  c  C  bulk  +  c  Si  surf  + 
cCsurf  +  cSibulk  +  c_SiH2_surf  +  c_CH2) 

Reactor  volume 

V_r 

Volumetric  flow  rate 

1 

Old  volumetric  flow 
rate 

1 

Volumetric  production 
rate 

R  const*T*  Vr*(r  1+r  2  +  r  3  +  r  4-r  7  -  r  8+r  9  +  r  10  +  r  11 
+  r_14)/p 

Old  volumetric 
production  rate 

R  const*T*  Vr*(r  1+r  2  +  r  3  +  r  4-r  7-r  8+r  9  +  r  10  +  r  11 
+  r_14)/p 

Reaction  sequence 
number 

17 

Sequence  number 

26 

Surface  reaction 
sequence  number 

15 

Temperature 

User  defined 

Temperature 

T_o 

Pressure 

User  defined 

Pressure 

R  const*T*(c  SiH4  +  c  SiH2  +  c  H2  +  c  Si2H6  +  c  Si  +  c  H  + 
c  C3H8  +  C  CH3  +  c  C2H5  +  c  CH4  +  c  C2H6  +  c  C2H4  + 
c_C2H2  +  c_Si2  +  c_Si2C  +  c_SiCH2  +  c_CH2) 

Importing... 

0 

Calculate 

thermodynamic 

properties 

On 

Calculate  transport 
properties 

On 

Reactor  type 

Batch,  constant  volume 

Mixture 

Gas 
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Description 

Value 

Solve  for  reactions  of 
type 

Volumetric 

Species  numbering 

Off 

Solvent  is  set 

Off 

Uniform  scaling  of 
concentration  variables 

0 

Show  equation 
assuming 

stdl/time 

2.2.1.  1:  SiH4<=>SiH2+H2 


The  settings  for  the  first  homogeneous  reaction  is  given.  Subsequent  reactions  follow  a  similar 
method  but  with  the  reaction-specific  data  listed  in  Appendix  B. 

Selection 


Geometric  entity  level  Entire  model 

Settings 

Description 

Value 

Formula 

SiH4<=>SiH2  +  H2 

Valid  reaction 

1 

Reaction  defined 

0 

Third  body 

Reset  to  default 

0 

Use  Arrhenius  expressions 

One 

Forward  frequency  factor 

1.67E-6 

Reverse  frequency  factor 

2.17E8 

Forward  temperature  exponent 

18.23 

Reverse  temperature  exponent 

3.1 

Forward  rate  constant 

Af  l*TAnf  l*exp(-Ef  1/(R  const*T)) 
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Description 

Value 

Reverse  rate  constant 

Ar  l*TAnr  l*exp(-Er  1/(R  const*T)) 

Equilibrium  constant 

1 

Specify  equilibrium  constant 

Zero 

Equilibrium  expression 

1.0 

Equilibrium  expression  numerator 

0 

Equilibrium  expression  denominator 

1.0 

Reaction  rate 

kf  l*c  SiE14-kr  l*c  SiEI2*c  E12 

Enthalpy  of  reaction 

-h_SiH4  +  h_SiH2  +  h_H2 

Entropy  of  reaction 

-s_SiH4  +  s_SiH2  +  s_H2 

Eleat  source  of  reaction 

-r_l*H_l 

Forward  stoichiometric  array 

{-1,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 

0,  0,  0} 

Reverse  stoichiometric  array 

{0,  1,  1,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 
0,0} 

Stoichiometric  array 

{-1,  1,  1,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 

0, 0,  0} 

Reactant  species 

SiH4 

Product  species 

{SiH2,  H2} 

Species 

{SiH4,  SiH2,  H2} 

Reaction  sequence  number 

1 

Reaction  type 

Reversible 

Forward  activation  energy 

-6.26E4 

Reverse  activation  energy 

-1.06E4 

Reaction  S 

0 

Reaction  Q 

0 

Reset  to  default 

0 

Reset  to  default 

0 
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Description 

Value 

Reset  to  default 

0 

Reset  to  default 

0 

Reset  to  default 

0 

Reaction  rate 

User  defined 

Enthalpy  of  reaction 

User  defined 

Entropy  of  reaction 

Automatic 

Heat  source  of  reaction 

Automatic 

Equilibrium  constant 

User  defined 

Rate 

0 

Equilibrium  Settings 

0 

Forward  rate  constant 

0 

Reverse  rate  constant 

0 
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2.2.2.  2:  Si2H6<=>SiH2+SiH4 

2.2.3.  3:  SiH2<=>Si+H2 

2.2.4.  5:  C3H8<=>CH3+C2H5 

2.2.5.  6:  CH4+H<=>CH3+H2 

2.2.6.  7:  C2H5+H<=>2CH3 

2.2.7.  8:  2CH3<=>C2H6 

2.2.8.  9:  C2H4+H<=>C2H5 

2.2.9.  10:  C2H4<=>C2H2+H2 

2.2.10.  13:  Si2<=>2Si 

2.2.11.  14:  Si2+CH4<=>Si2C+2H2 

2.2.12.  15:  SiH2+Si<=>Si2+H2 

2.2.13.  16:  CH3+Si<=>SiCH2+H 

2.2.14.  17:  SiCH2+SiH2<=>Si2C+2H2 


2.2.15.  1:  C(b)+Si(s)+H2=>SiH2+C(s) 


The  settings  for  the  first  surface  reaction  is  given.  Subsequent  surface  reactions  follow  a  similar 
method  but  with  the  reaction-specific  data  taken  from  Blanquet,  et  al. 

Selection 


Geometric  entity  level 


Entire  model 


Settings 


Description 

Value 

Formula 

C(b)  +  Si(s)  +  H2=>SiH2  +  C(s) 

Valid  reaction 

1 

Reaction  defined 

0 

Reset  to  default 

0 

Reaction  type 

Irreversible 

Third  body 

140 


Chemical  Reaction  Model 


Description 

Value 

Use  Arrhenius  expressions 

1 

Forward  sticking  coefficient 
frequency  factor 

2.200el8 

Forward  sticking  coefficient 
exponent 

0 

Forward  sticking  coefficient 
activation  energy 

52786[K]*R 

Forward  rate  constant 

af  l*TAbf  l*exp(-ef  1/(R  const*T)) 

Reaction  rate 

ksf  l*c  FI2*c  C  bulk*c  Si  surf 

Enthalpy  of  reaction 

h_SiH2  -  h_H2  -  h  C  bulk  -  h  Si  surf  +  h  C  surf 

Entropy  of  reaction 

s_SiFI2  -  s_H2  -  s  C  bulk  -  s  Si  surf  +  s  C  surf 

Surface  heat  source  of  reaction 

-q _ 1  *H  1 

Forward  stoichiometric  array 

{0,  0,  -1,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  -1,  -1,  0,  0, 
0,0} 

Reverse  stoichiometric  array 

{0,  1,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  1,  0,  0, 
0} 

Stoichiometric  array 

{0,  1,  -1,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  -1,  -1,  1,  0, 
0,0} 

Reactant  species 

(C(b),  Si(s),  H2} 

Product  species 

{SiH2,  C(s)} 

Species 

{C  bulk,  Si  surf,  F12,  SiF12,  C  surf} 

Reaction  sequence  number 

1 
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2.2.16.  2:  2Si(b)+2C(s)+H2=>C2H2+2Si(s) 

2.2.17.  3:  SiH4+C(s)=>SiH2(s)+H2+C(b) 

2.2.18.  4:  SiH2(s)=>H2+Si(s) 

2.2.19.  5:  SiH2+C(s)=>SiH2(s)+C(b) 

2.2.20.  6:  Si+C(s)=>Si(s)+C(b) 

2.2.21.  7:  C2H2+Si(s)=>2C(s)+H2+2Si(b) 

2.2.22.  8:  C2H4+2Si(s)=>2C(s)+2H2+2Si(b) 

2.2.23.  9:  CH4+Si(s)=>C(s)+2H2+Si(b) 

2.2.24.  11:  CH3+Si(s)=>C(s)+1.5H2+Si(b) 

2.2.25.  12:  Si2+2C(s)=>2Si(s)+2C(b) 

2.2.26.  13:  Si2C+Si(s)=>Si2+C(s)+Si(b) 

2.2.27.  14:  SiCH2+C(s)=>Si(s)+CH2+C(b) 

2.2.28.  15:  CH2+Si(s)=>C(s)+H2+Si(b) 

2.2.29.  Species:  SiH4 

The  settings  for  the  first  species  is  given.  Subsequent  species  follow  a  similar  method  but  with  the 
species-specific  data  listed  in  Appendix  A. 

Selection 

Geometric  entity  level  Entire  model 


Settings 

Description 

Value 

Species  is  defined 

1 

Dependent 

0 

Sequence  number 

2 

Lock  concentration/ 
activity 

0 

Species  name 

SiH4 
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Description 

Value 

Species  label 

SiH4 

Reset  to  default 

0 

Molecular  weight 

0.0321 17[kg/mol] 

Charge 

0 

Old  rate  expression 

-r_l  +  r_2 

Rate  expression 

-r  1  +  r_2 

Surface  rate  expression 

-q_3 

Old  surface  rate 
expression 

-q_3 

Lower  temperature  limit 

298  [K] 

Middle  temperature  limit 

1300[K] 

Upper  temperature  limit 

6000[K] 

Polynomial  coefficients 

{6.060189,  139.9632,  -77.88474,  16.24095,  0.135509,  -6.91824, 
174.3351} 

Polynomial  coefficients 

{99.84949,  4.25153,  -0.809269,  0.053437,  -20.39005,  -74.8492, 
266.8015} 
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Description 

Molar  entropy 


Molar  enthalpy 


Value 


(T<=Tlo_SiH4)*l[J/ 

( mol*K)]  *(  aLo  1  _SiH4*log(Tlo_SiH4/ 1 000  [K] )  + 
aLo2_SiH4*(Tlo_SiH4/l  000)  + 

( l/2)*aLo3_S  iH4  *  (Tlo_S  iH4/ 1 000)A2  + 

( 1  /3 ) * aLo4_S iH4 * (Tlo_S iH4/ 1 000)A3  - 

(l/2)*aLo5_SiH4*l[KA6]*(Tlo_SiH4/1000)A  -  2  +  aLo7_SiH4)  + 
(T>Tlo_SiH4)*(T<=Tmid_SiH4)*  1  [J/ 

(mol*K)]*(aLol_SiH4*log(T/1000[K])  +  aLo2_SiH4*(T/1000)  + 
(l/2)*aLo3_SiH4*(T/1000)A2  +  (l/3)*aLo4_SiH4*(T/1000)A3  - 
( l/2)*aLo5_SiH4*l[KA6]*(T/1000)A  -  2  +  aLo7_SiH4)  + 

(T>T  mid_SiH4)  *  (T  <=Thi_SiH4)  *  1  [J/ 

(mol*K)]*(aHil_SiH4*log(T/1000[K])  +  aHi2_SiH4*(T/1000)  + 

( l/2)*aHi3_SiH4*(T/l 000)A2  +  (l/3)*aHi4_SiH4*(T/1000)A3  - 
( l/2)*aHi5_SiH4*  1  [KA6]*(T/1 000)A  -  2  +  aHi7_SiH4)  + 
(T>Thi_SiH4)*l[J/(mol*K)]*(aHil_SiH4*log(Thi_SiH4/1000[K]) 
+  aHi2_SiH4*(Thi_SiH4/1000)  + 

( l/2)*aHi3_S  iH4  *(Thi_S  iH4/ 1 000)A2  + 

( 1/3)* aHi4_S iH4 *(Thi_S iH4/ 1 000)A3  - 
( !/2)*aHi5_SiH4*  1  [KA6]*(Thi_SiH4/l  000)A  -  2  +  aHi7_SiH4) 


(T<=Tlo_SiH4)*l[J/(mol*K)]*(aLol_SiH4*(Tlo_SiH4/1000)  + 

( l/2)*aLo2_SiH4*((Tlo_SiH4/1000)A2)  + 

( l/3)*aLo3_SiH4*((Tlo_SiH4/1000)A3)  + 

( l/4)*aLo4_SiH4*((Tlo_SiH4/1000)A4)  - 

aLo5_SiH4*  1  [KA6]*((Tlo_SiH4/l 000)A  -  1)  +  aLo6_SiH4)  + 

(T>Tlo_SiH4)*(T<=Tmid_SiH4)*l[J/(mol*K)]*(aLol_SiH4*(T/ 

1000)  +  (l/2)*aLo2_SiH4*((T/1000)A2)  +  (l/3)*aLo3_SiH4*((T/ 

1000)A3)  +  (l/4)*aLo4_SiH4*((T/1000)A4)  - 

aLo5_SiH4*  1  [KA6]*((T/1 000)A  -  1)  +  aLo6_SiH4)  + 

(T>T  mid_SiH4)  *  (T  <=Thi_SiH4)  *  1  [J/(mol*K)]*(aHil_SiH4*(T/ 
1000)  +  (l/2)*aHi2_SiH4*((T/1000)A2)  +  (l/3)*aHi3_SiH4*((T/ 
1000)A3)  +  (l/4)*aHi4_SiH4*((T/1000)A4)  - 
aHi5_SiH4*l[KA6]*((T/1000)A  -  1)  +  aHi6_SiH4)  + 
(T>Thi_SiH4)*l[J/(mol*K)]*(aHil_SiH4*(Thi_SiH4/1000)  + 
(l/2)*aHi2_SiH4*((Thi_SiH4/1000)A2)  + 
(l/3)*aHi3_SiH4*((Thi_SiH4/1000)A3)  + 
(l/4)*aHi4_SiH4*((Thi_SiH4/1000)A4)  - 
aHi5_SiH4*  1  [KA6]*((Thi_SiH4/l  000)A  -  1)  +  aHi6_SiH4) 
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Description  Value 

Heat  capacity  at  constant 
pressure 

(T<=Tlo  SiH4)*l[J/(mol*K)]*(aLol  SiH4  + 

aLo2  SiH4*(Tlo  SiH4/1000)  +  aLo3  SiH4*(Tlo  SiH4/1000)A2  + 

aLo4  SiH4*(Tlo  SiH4/1000)A3  + 

aLo5  SiH4*  1  [KA6]*(Tlo  SiH4/1000)A  -  2)  + 

(T>Tlo  SiH4)*(T<=Tmid  SiH4)*l[J/(mol*K)]*(aLol  SiH4 + 
aLo2  SiH4*(T/1000)  +  aLo3  SiH4*(T/1000)A2  +  aLo4  SiH4*(T/ 
1000)A3  +  aLo5  SiH4*l[KA6]*(T/1000)A  -  2)  + 

(T>Tmid  SiH4)*(T<=Thi  SiH4)*l[J/(mol*K)]*(aHil  SiH4  + 
aHi2  SiH4*(T/1000)  +  aHi3  SiH4*(T/1000)A2  +  aHi4  SiH4*(T/ 
1000)A3  +  aHi5  SiH4*l[KA6]*(T/1000)A  -  2)  +  (T>Thi  SiH4)*l[J/ 
(mol*K)]*(aHil  SiH4  +  aHi2  SiH4*(Thi  SiH4/1000)  + 
aHi3  SiH4*(Thi  SiH4/1000)A2  +  aHi4  SiH4*(Thi  SiH4/1000)A3 
+  aHi5_SiH4*  1  [KA6] *(Thi_SiH4/l  000)A  -  2) 

Potential  characteristic 
length 

3.458[angstrom] 

Potential  energy 
minimum 

107.4[K] 

Dipole  moment 

0 

Diffusivity 

2.6950000000000004E-22*sqrt(5.0E-4*(M_SiH4  +  M_H2)*TA3/ 

(M  SiH4*M  H2))/(p*sigma  SiH4*sigma  H2*(  1.06036/(T/ 
sqrt(epsilonkb  SiH4*epsilonkb  H2))A0.1561  +  0.193*exp(- 
(0.47635*T/sqrt(epsilonkb_SiH4*epsilonkb_H2)))  +  1.03587*exp(- 
(1.52996*T/sqrt(epsilonkb  SiH4*epsilonkb  H2)))  +  1.76474*exp(- 
(3.8941 1  *T/sqrt(cpsilonkb_SiI  I4*epsilonkb_1 12))))) 

Mobility 

re.D  SiH4*e  const/k  B  const/T 

Dynamic  viscosity 

2.669e-6*sqrt(T*M  SiH4*le3)/((sigma  SiH4)*lelO)A2*(1.16145/ 
(T/epsilonkb_SiH4)A0. 14874  +  0.52487/exp(0.77320*T/ 
epsilonkb  SiH4)  +  2.16178/exp(2.43787*T/epsilonkb  SiH4)  + 
4.998e-40*mu_SiH4A4/ 

(((sigma  SiH4)*lelO)A6*k  B  constA2*epsilonkb  SiH4*T))A  -  1 

Thermal  conductivity 

eta  SiH4/M  SiH4*(1.15*Cp  SiH4  +  0.88*R  const) 

145 


Chemical  Reaction  Model 


Description 

Value 

Feed  stream  molar 
enthalpy 

R  const*((Tf  SiFI4<=Tlo  SiF14)*Tlo  SiF14*(aLol  SiF14  + 

0.5*aLo2  SiH4*Tlo  SiH4  +  aLo3  SiH4*Tlo  SiH4A2/3  + 

0.25*aLo4  SiH4*Tlo  SiH4A3  +  0.2*aLo5  SiH4*Tlo  SiH4A4  + 
aLo6  SiH4/Tlo  SiH4)  + 

(Tf  SiH4>Tlo  SiH4)*(Tf  SiH4<=Tmid  SiH4)*Tf  SiH4*(aLol  Si 
H4  +  0.5*aLo2  SiH4*Tf  SiH4  +  aLo3  SiH4*Tf  SiH4A2/3  + 
0.25*aLo4  SiH4*Tf  SiH4A3  +  0.2*aLo5  SiH4*Tf  SiH4A4  + 
aLo6  SiH4/Tf  SiH4)  + 

(Tf  SiH4>Tmid  SiH4)*(Tf  SiH4<=Thi  SiII4)*Tf  SiH4*(aHil  Si 
H4  +  0.5*aHi2  SiH4*Tf  SiH4  +  aHi3  SiH4*Tf  SiH4A2/3  + 
0.25*aHi4  SiH4*Tf  SiH4A3  +  0.2*aHi5  SiH4*Tf  SiH4A4  + 
aHi6  SiH4/Tf  SiH4)  + 

(Tf  SiH4>Thi  SiH4)*Thi  SiH4*(aHil  SiH4  + 

0.5*aHi2  SiH4*Thi  SiH4  +  aHi3  SiH4*Thi  SiH4A2/3  + 

0.25*aHi4  SiH4*Thi  SiH4A3  +  0.2*aHi5  SiH4*Thi  SiH4A4  + 
aHi6_SiH4/Thi_SiH4)) 

Volumetric  feed  rate 

0 

Feed  stream  temperature 

0 

Additional  enthalpy 
contribution 

0 

Species  type 

None 

Density 

c_SiH4*M_SiH4 

Initial  concentration 

c_SiF14_o 

Specify  diffusivity  and 
mobility 

1 

Rate  expression 

Automatic 

Surface  rate  expression 

User  defined 

Species  rate 

0 

Reset  to  default 

0 

Reset  to  default 

0 
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2.2.30. 

Species:  SiH2 

2.2.31. 

Species:  H2 

2.2.32. 

Species:  Si2H6 

2.2.33. 

Species:  Si 

2.2.34. 

Species:  H 

2.2.35. 

Species:  C3H8 

2.2.36. 

Species:  CH3 

2.2.37. 

Species:  C2H5 

2.2.38. 

Species:  CH4 

2.2.39. 

Species:  C2H6 

2.2.40. 

Species:  C2H4 

2.2.41. 

Species:  C2H2 

2.2.42. 

Species:  Si2 

2.2.43. 

Species:  Si2C 

2.2.44. 

Species:  SiCH2 

2.2.45. 

Species:  C(b) 

2.2.46. 

Species:  Si(s) 

2.2.47. 

Species:  C(s) 

2.2.48. 

Species:  Si(b) 

2.2.49. 

Species:  SiH2(s) 

2.2.50. 

Species:  CH2 
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3.  Study  1 

3.1.  Parametric  Sweep 

Parameter  name:  T  o 

Parameters: 

3.2.  Time  Dependent 
Study  settings 


Property 

Value 

Include  geometric  nonlinearity 

Off 

Times:  range(0,le-8,2e-6) 

Physics  selection 


Physics 

Discretization 

Reaction  Engineering  (re) 

physics 

3.3.  Solver  Configurations 
3.3.1.  Solver  1 

Compile  Equations:  Time  Dependent  (stl) 
Study  and  step 


Name 

Value 

Use  study 

Study  1 

Use  study  step 

Time  Dependent 

Dependent  Variables  1  (vl) 
General 


Name 

Value 

Defined  by  study  step 

Time  Dependent 

Initial  values  of  variables  solved  for 


Name 

Value 

Solution 

Zero 
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Values  of  variables  not  solved  for 


Name 

Value 

Solution 

Zero 

Concentration  (modl.ODEl)  (modlODEl) 

General 


Name 

Value 

State  components 

modl.re.c  SiH4 

Concentration  (modl.ODE2)  (modl_ODE2) 

General 


Name 

Value 

State  components 

modl.re.c  SiH2 

Concentration  (modl.ODE3)  (modl_ODE3) 

General 


Name 

Value 

State  components 

modl.re.c  H2 

Concentration  (modl.ODE4)  (modl_ODE4) 

General 


Name 

Value 

State  components 

modl.re.c_Si2H6 

Concentration  (modl.ODE5)  (modi  ODE5) 

General 

Name 

Value 

State  components 

modl.re.c_Si 

Concentration  (modl.ODE6)  (modl_ODE6) 

General 
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Name 

Value 

State  components 

modi. re. c  H 

Concentration  (modl.ODE7)  (modl_ODE7) 

General 


Name 

Value 

State  components 

modl.re.c_C3H8 

Concentration  (modl.ODE8)  (modl_ODE8) 

General 


Name 

Value 

State  components 

modl.re.c_CH3 

Concentration  (modl.ODE9)  (modl_ODE9) 

General 


Name 

Value 

State  components 

modi. re. c  C2H5 

Concentration  (modl.ODElO)  (modl  ODElO) 

General 


Name 

Value 

State  components 

modi. re. c  CH4 

Concentration  (modl.ODEll)  (modl  ODEll) 

General 


Name 

Value 

State  components 

modi. re. c  C2H6 

Concentration  (modl.ODE12)  (modl_ODE12) 

General 
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Name 

Value 

State  components 

modl.re.c  C2H4 

Concentration  (modl.ODE13)  (modl_ODE13) 

General 


Name 

Value 

State  components 

modl.re.c  C2H2 

Concentration  (modl.ODE14)  (modl_ODE14) 

General 


Name 

Value 

State  components 

modl.re.c_Si2 

Concentration  (modl.ODE15)  (modl_ODE15) 

General 


Name 

Value 

State  components 

modl.re.c_Si2C 

Concentration  (mod l.ODE  16)  (modl_ODE16) 

General 


Name 

Value 

State  components 

modl.re.c_SiCH2 

Concentration  (modl.ODE17)  (modl_ODE17) 

General 


Name 

Value 

State  components 

mod  1  .re.cCbulk 

Concentration  (modl.ODE18)  (modl_ODE18) 

General 
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Name 

Value 

State  components 

modi  .re.cSisurf 

Concentration  (modl.ODE19)  (modl_ODE19) 

General 


Name 

Value 

State  components 

mod  1  .re.cCsurf 

Concentration  (modl.ODE20)  (modl_ODE20) 

General 


Name 

Value 

State  components 

modi  .re.cSibulk 

Concentration  (modl.ODE21)  (modl_ODE21) 

General 


Name 

Value 

State  components 

modl.re.c  CH2 

Time-Dependent  Solver  1  (tl) 

General 

Name 

Value 

Defined  by  study 
step 

User  defined 
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Name 

Value 

Time 

{0,  1.0E-8,  2.0E-8,  3.0000000000000004E-8,  4.0E-8,  5.0E-8, 
6.00000000000000 IE-8,  7.0E-8,  8.0E-8,  9.0E-8,  1.0E-7,  1.1E-7, 
1.2000000000000002E-7,  1.3E-7,  1.4E-7,  1.5E-7,  1.6E-7, 

1. 700000000000000 IE-7,  1.8E-7,  1.9E-7,  2.0E-7,  2.1E-7,  2.2E-7,  2.3E-7, 
2.4000000000000003E-7,  2.5E-7,  2.6E-7,  2.7E-7,  2.8E-7, 
2.9000000000000003E-7,  3.0E-7,  3.1E-7,  3.2E-7,  3.3E-7, 
3.4000000000000003E-7,  3.5E-7,  3.6E-7,  3.7E-7,  3.8E-7,  3.9E-7,  4.0E-7, 
4.1E-7,  4.2E-7,  4.3E-7,  4.4E-7,  4.5000000000000003E-7,  4.6E-7,  4.7E-7, 

4. 80000000000000 IE-7,  4.9E-7,  5.0E-7,  5.1E-7,  5.2E-7,  5.3E-7,  5.4E-7, 
5.5E-7,  5.6E-7,  5.7E-7,  5.800000000000001E-7,  5.9E-7,  6.0E-7,  6.1E-7, 
6.2E-7,  6.3E-7,  6.4E-7,  6.5E-7,  6.6E-7,  6.7E-7,  6.800000000000001E-7, 
6.90000000000000 IE-7,  7.0E-7,  7.1E-7,  7.2E-7,  7.3E-7,  7.4E-7,  7.5E-7, 
7.6E-7,  7.7E-7,  7.8E-7,  7.900000000000001E-7,  8.0E-7,  8.1E-7,  8.2E-7, 
8.3E-7,  8.4E-7,  8.5E-7,  8.6E-7,  8.7E-7,  8.8E-7,  8.900000000000001E-7, 
9.00000000000000 IE-7,  9.1E-7,  9.2E-7,  9.3E-7,  9.4E-7,  9.5E-7, 
9.60000000000000 IE-7,  9.7E-7,  9.8E-7,  9.9E-7,  1.0E-6} 

Relative  tolerance 

le-9 

Absolute  tolerance 


Name 

Value 

Tolerance 

1.0E-4 

Output 


Name 

Value 

Times  to  store 

Steps  taken  by  solver 

Advanced 


Name 

Value 

Consistent  initialization 

Off 

Name 

Value 

Solution 

Solver  1 
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1.  Global  Definitions 
1.1.  Parameters  1 
Parameters 


Name 

Expression 

Description 

R  fib 

0.005[m] 

fiber  radius 
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2.  Model  1  (modi) 

2.1.  Definitions 

2.1.1.  Variables 
Variables  la 
Selection 

Geometric  entity  level  Entire  model 


Name 

Expression 

Description 

density 

10[mol/mA3] 

material  density 

vx  int 

-(chds  ,ntflux_c  *  chds  .nx)/  density 

surface  velocity  in  x-direction 

vy  int 

-(chds.ntflux_c*chds.ny)/density 

surface  velocity  in  y-direction 

2.1.2.  Coordinate  Systems 
Boundary  System  1 


Coordinate  system  type 

Boundary  system 

Identifier 

sysl 

Settings 


Name 

Value 

Coordinate  names 

{tl,  n,  to} 

Create  first  tangent  direction  from 

Global  Cartesian  (spatial) 
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2.2.  Geometry  1 


Geometry  1 

Units 


Length  unit 

m 

Angular  unit 

deg 

Geometry  statistics 


Property 

Value 

Space  dimension 

2 

Number  of  domains 

1 

Number  of  boundaries 

20 

2.2.1.  Square  1  (sql) 
Position 


Name 

Value 

Position 

{0,0} 
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Name 

Value 

Side  length 

0.1 

Side  length 

0.1 

2.2.2.  Circle  1  (cl) 
Position 


Name 

Value 

Position 

{0.05,  0.05} 

Radius 

R  fib 

2.2.3.  Circle  2  (c2) 
Position 


Name 

Value 

Position 

{0.05,  0.07} 

Radius 

R_fib 

2.2.4.  Circle  3  (c3) 
Position 


Name 

Value 

Position 

{0.02,  0.03} 

Radius 

R  fib 

2.2.5.  Circle  4  (c4) 
Position 


Name 

Value 

Position 

{0.06,  0.04} 

Radius 

R  fib 

2.2.6.  Compose  1  (col) 
Selections  of  resulting  entities 


Moving  Mesh  Model 


Name 

Value 

Create  selections 

On 

Set  formula 

sql  -  cl  -  c2  -  c3  -  c4 

2.3.  Transport  of  Diluted  Species  (chds) 

or 

009 

out 

O.OT 
0  06 
OH 

oof 
0  01 
0  07 

oor 

0 

•0  02  0  0  02  0  04  0  0ft  00ft  01  0 


Transport  of  Diluted  Species 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Concentration 

Linear 

Compute  boundary  fluxes 

1 

Apply  smoothing  to  boundary  fluxes 

1 
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Description 

Value 

Value  type  when  using  splitting  of  complex 
variables 

Real 

Equation  form 

Study  controlled 

Migration  in  electric  field 

0 

Convection 

0 

Convective  term 

Non  -  conservative  form 

Equation  residual 

Approximate  residual 

Streamline  diffusion 

1 

Crosswind  diffusion 

1 

Crosswind  diffusion  type 

Do  Carmo  and  Galeao 

Enable  space-dependent  physics  interfaces 

0 

Synchronize  with  COMSOL  Multiphysics 

Show  equation  assuming 

stdl/time 
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Diffusion 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Velocity  field 

User  defined 

Velocity  field 

{0,  0,  0} 

Electric  potential 

User  defined 

Electric  potential 

0 

Diffusion  coefficient 

User  defined 

Diffusion  coefficient 

{{le-9[mA2/s],  0,  0},  {0,  le-9[mA2/s],  0},  {0,  0,  le-9[mA2/s]}} 

Bulk  material 

None 
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2.3.2.  No  Flux  1 


No  Flux  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundaries  1-4 

Settings 


Description 

Value 

Apply  for  all  species 

Apply  for  all  species 
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Initial  Values  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Concentration 

20 
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2.3.4.  Concentration  1 


0  01  0  02  0  0)  0  04  0  0)  0  0«  0  07 


Concentration  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundaries  5-20 

Settings 


Description 

Value 

Concentration 

0 

Species  c 

1 

Apply  reaction  terms  on 

All  physics  (symmetric) 

Use  weak  constraints 

0 
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2.4.  Deformed  Geometry  (dg) 

or 

009 
OOt 

o.o7 
006 
0  05 

oo<r 

0  09 
0  0? 
oor 
0 

■0  02  0  0  02  0  04  0  06  0  06  0  1  0 


Deformed  Geometry’ 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Equation  form 

Study  controlled 

Geometry  frame  coordinates 

{Xg,  Yg,  Zg} 

Geometry  shape  order 

1 

Mesh  smoothing  type 

Laplace 

Show  equation  assuming 

0 

2.4.1.  Fixed  Mesh  1 
Selection 
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Geometric  entity  level 

Domain 

Selection 

No  domains 

2.4.2.  Prescribed  Mesh  Displacement  1 


Prescribed  Mesh  Displacement  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundaries  2-3 

Settings 


Description 

Value 

Prescribed  #  displacement 

{0,1} 

Prescribed  mesh  displacement 

{0,0} 

Use  weak  constraints 

0 
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Free  Deformation  1 

Selection 


Geometric  entity  level 

Domain 

Selection 

Domain  1 

Settings 


Description 

Value 

Initial  mesh  displacement 

{0,0} 
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2.4.4.  Prescribed  Mesh  Displacement  2 


or 

009 

0  M 

o  oT 

0  Ot 

o 

0  OS 

o 

0  04 

■  o 

0  03 

0  02 

o 

oof 

0 


■0  02  0  0  02  0  04  0  00 


Prescribed  Mesh  Displacement  2 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundaries  1,  4 

Settings 


Description 

Value 

Prescribed  #  displacement 

{1,0} 

Prescribed  mesh  displacement 

{0,0} 

Use  weak  constraints 

0 

0  00  0  1  0 
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2.4.5.  Prescribed  Mesh  Velocity  1 


C  075" 
0  07 
0  0»>S 

0  u* 


C  051 
0  04 
e  04b 
004 
C  034 

OOJ 
C  025 


001  002 


003 


0  04  0  04  0  04 


Prescribed  Mesh  Velocity  1 

Selection 


Geometric  entity  level 

Boundary 

Selection 

Boundaries  5-20 

Settings 


Description 

Value 

Prescribed  #  velocity 

{1,1} 

Prescribed  mesh  velocity 

{vx_int,  vy_int} 

Use  weak  constraints 

0 

2.5.  Meshes 
2.5.1.  Mesh  1 
Mesh  statistics 


0  07 
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Property 

Value 

Minimum  element  quality 

0.6886 

Average  element  quality 

0.9202 

Triangular  elements 

2364 

Edge  elements 

172 

Vertex  elements 

20 

or 

0  09 
0  M 

0.07 
006 
0  06 
0  04" 

0  01 
0  0? 

OOf 
0 

0  02 


0  02 


004 


006 


006 


Mesh  1 

Size  (size) 

Settings 


Name 

Value 

Maximum  element  size 

0.0053 

Minimum  element  size 

3.0E-5 

Resolution  of  curvature 

0.3 

Maximum  element  growth  rate 

1.3 

Predefined  size 

Fine 
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Free  Triangular  1  (ftril) 
Selection 

Geometric  entity  level 


Remaining 
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3.  Study  1 

3.1.  Time  Dependent 

Study  settings 


Property 

Value 

Include  geometric  nonlinearity 

Off 

Times:  range(0,100,10A6) 

Mesh  selection 


Geometry 

Mesh 

Geometry  1  (geoml) 

meshl 

Physics  selection 


Physics 

Discretization 

Transport  of  Diluted  Species  (chds) 

physics 

Deformed  Geometry  (dg) 

physics 

_ 

3.2.  Solver  Configurations 
3.2.1.  Solver  1 

Compile  Equations:  Time  Dependent  (stl) 
Study  and  step 


Name 

Value 

Use  study 

Study  1 

Use  study  step 

Time  Dependent 

Dependent  Variables  1  (vl) 
General 


Name 

Value 

Defined  by  study  step 

Time  Dependent 
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Initial  values  of  variables  solved  for 


Name 

Value 

Solution 

Zero 

Values  of  variables  not  solved  for 


Name 

Value 

Solution 

Zero 

Concentration  (modl.c)  (modl_c) 


General 

Name 

Value 

Field  components 

modl.c 

modl.xy  (modl_xy) 

General 


Name 

Value 

Field  components 

{x,y} 

modl.XgYg  (modl_XgYg) 

General 

Name 

Value 

Field  components 

{Xg,  Yg} 

Solve  for  this  field 

Off 

Time-Dependent  Solver  1  (tl) 
General 


Name 

Value 

Defined  by  study 
step 

Time  Dependent 
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Name 

Time 


Value 


{0,  100,  200,  300,  400,  500,  600,  700,  800,  900,  1000,  1100,  1200,  1300, 
1400,  1500,  1600,  1700,  1800,  1900,  2000,  2100,  2200,  2300,  2400,  2500, 

2600,  2700,  2800,  2900,  3000,  3100,  3200,  3300,  3400,  3500,  3600,  3700, 

3800,  3900,  4000,  4100,  4200,  4300,  4400,  4500,  4600,  4700,  4800,  4900, 

5000,  5100,  5200,  5300,  5400,  5500,  5600,  5700,  5800,  5900,  6000,  6100, 

6200,  6300,  6400,  6500,  6600,  6700,  6800,  6900,  7000,  7100,  7200,  7300, 

7400,  7500,  7600,  7700,  7800,  7900,  8000,  8100,  8200,  8300,  8400,  8500, 

8600,  8700,  8800,  8900,  9000,  9100,  9200,  9300,  9400,  9500,  9600,  9700, 

9800,  9900,  10000,  10100,  10200,  10300,  10400,  10500,  10600,  10700, 


10800,  10900, 11000,  11100,  11200, 
11800,  11900,  12000,  12100,  12200, 
12800,  12900,  13000,  13100,  13200, 
13800, 13900, 14000,  14100,  14200, 
14800, 14900, 15000,  15100, 15200, 
15800, 15900, 16000, 16100, 16200, 
16800, 16900, 17000, 17100, 17200, 
17800,  17900, 18000,  18100,  18200, 
18800,  18900, 19000, 19100, 19200, 
19800,  19900, 20000, 20100, 20200, 
20800,  20900, 21000,  21100,  21200, 
21800,  21900, 22000,  22100,  22200, 
22800,  22900,  23000,  23100,  23200, 
23800, 23900, 24000,  24100, 24200, 
24800, 24900, 25000, 25100, 25200, 
25800,  25900, 26000,  26100,  26200, 
26800,  26900, 27000,  27100, 27200, 
27800,  27900, 28000,  28100,  28200, 
28800,  28900,  29000,  29100,  29200, 
29800,  29900, 30000,  30100,  30200, 
30800,  30900,  31000,  31100,  31200, 
31800, 31900, 32000, 32100, 32200, 
32800, 32900, 33000, 33100, 33200, 
33800,  33900, 34000,  34100,  34200, 
34800,  34900, 35000,  35100,  35200, 
35800,  35900, 36000,  36100, 36200, 
36800,  36900, 37000,  37100, 37200, 
37800,  37900, 38000,  38100,  38200, 
38800,  38900,  39000,  39100,  39200, 
39800,  39900, 40000,  40100,  40200, 
40800,  40900,  41000,  41100,  41200, 
41800, 41900, 42000,  42100, 42200, 
42800, 42900, 43000, 43100, 43200, 
43800,  43900, 44000,  44100,  44200, 

44800,  4490Q, 45000.  45100. 45200. 
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11300,  11400,  11500,  11600,  11700, 
12300,  12400,  12500,  12600,  12700, 
13300, 13400,  13500,  13600,  13700, 
14300, 14400,  14500,  14600,  14700, 
15300, 15400,  15500,  15600,  15700, 
16300, 16400,  16500,  16600,  16700, 
17300, 17400,  17500,  17600,  17700, 
18300, 18400,  18500,  18600,  18700, 
19300, 19400,  19500,  19600,  19700, 
20300, 20400,  20500,  20600,  20700, 
21300,  21400,  21500,  21600,  21700, 
22300, 22400,  22500,  22600,  22700, 
23300, 23400,  23500,  23600,  23700, 
24300, 24400,  24500,  24600,  24700, 
25300, 25400,  25500,  25600,  25700, 
26300, 26400,  26500,  26600,  26700, 
27300, 27400,  27500,  27600,  27700, 
28300, 28400,  28500,  28600,  28700, 
29300, 29400,  29500,  29600,  29700, 
30300, 30400,  30500,  30600,  30700, 
31300,  31400,  31500,  31600,  31700, 
32300, 32400,  32500,  32600,  32700, 
33300, 33400,  33500,  33600,  33700, 
34300, 34400,  34500,  34600,  34700, 
35300, 35400,  35500,  35600,  35700, 
36300, 36400,  36500,  36600,  36700, 
37300, 37400,  37500,  37600,  37700, 
38300,  38400,  38500,  38600,  38700, 
39300, 39400,  39500,  39600,  39700, 
40300, 40400,  40500,  40600,  40700, 
41300,  41400,  41500,  41600,  41700, 
42300, 42400,  42500,  42600,  42700, 
43300, 43400,  43500,  43600,  43700, 
44300, 44400,  44500,  44600,  44700, 
45300.  45400.  45500.  45600.  45700. 
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General 


Name 

Value 

Remesh  in  geometry 

Geometry  1 

Output 


Name 

Value 

Solution 

Automatic  Remeshing  2 

Meshes 

{mesh2,  mesh3,  mesh4, 
mesh5,  mesh6} 

Fully  Coupled  1  (fcl) 

General 


Name 

Value 

Linear  solver 

Direct  1 

Automatic  Remeshing  2 
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4.  Results 

4.1.  Data  Sets 

4.1.1.  Solution  1 
Selection 


Geometric  entity  level 

Domain 

Selection 

Geometry  geoml 

Solution 


Name 

Value 

Solution 

Solver  1 

Model 

Save  Point  Geometry  1 

4.1.2.  Solution  2 
Selection 


Geometric  entity  level 

Domain 

Selection 

Geometry  geoml 

Solution 


Name 

Value 

Solution 

Automatic  Remeshing  2 

Model 

Save  Point  Geometry  1 
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4.2.  Plot  Groups 

4.2.1.  Concentration  (chds) 

(.OA«n*r*tiOA  (moVm1)  C6<K«*rr«T>on  (moi/vnh 


A  20  412 


0  0  02  0  04  0  06  0  06  01  WO 


Surface:  Concentration  (mol/m3)  Surface:  Concentration  (mol/m3) 
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